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SUMMARY 


PROBLEM 

Discuss,  in  a  form  suitable  for  comparison,  the  most  important  characteristics  of 
underwater  information-transfer  methods.  Tills  information  should  help  designers  inter¬ 
ested  in  underwater  communications,  navigation,  detection,  and  identification. 

RESULTS 


1 .  Tradeoffs  involved  in  the  undersea  detection  of  static  and  time-varying  signals 
were  examined,  and  quantitative  results  were  obtained. 

2.  Static  electric  and  magnetic  fields,  radioactive  fields,  chemical  systems, 
cxtremcly-low-frcqucncy  (ELF)  electromagnetic  radiation,  light,  and  acoustics  were 
studied,  and  the  results  were  summarized  in  chart  form. 

3.  For  time-varying  fields,  underwater  light  exhibits  inferior  range  and  ELF  exhibits 
inferior  bandwidth  when  compared  with  acoustics. 

4.  When  used  as  static  markers,  ranges  to  approximately  3  and  14  km  can  be  achieved 
for  chemical  and  ELF  sources,  respectively. 

5.  Considerations  beyond  the  preliminary'  design  stage  are  not  discussed  in  this 
report,  and  more  detailed  treatises,  many  of  which  are  referenced,  should  be  consulted. 
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A  —  ampere 

km 

A  —  angstrom 
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INTRODUCTION 


To  those  who  are  accustomed  to  air  as  the  supporting  or  propagating  medium,  the 
underwater  transfer  of  information  presents  new  problems:  In  the  ocean  environment, 
some  standard  information-transferring  techniques  which  arc  used  in  the  atmosphere  are 
nearly  useless,  and  other  techniques  which  arc  inappropriate  in  a  gaseous  environment  per¬ 
form  well.  For  example,  high-frequency  electromagnetic  radiation,  the  most  widely  used 
information  carrier  on  iand.  is  severely  attenuated  in  conducting  seawater:  acoustic  propa¬ 
gation.  however,  exhibits  lower  attenuation  in  the  sea  than  in  the  air. 

The  purpose  of  this  report  is  to  discuss,  in  a  form  suitable  for  comparison,  the 
various  techniques  that  arc  useful  in  the  underwater  transfer  of  information.  As  informa¬ 
tion  transfer  is  fundamental  to  problems  of  underwater  markers,  communication,  naviga¬ 
tion.  detection,  and  identification,  it  is  hoped  that  the  comparisons  in  th*s  report  will  be 
useful  to  a  wide  variety  of  investigators  who  arc  unfamiliar  with  the  special  limitations 
of  the  ocean  environment.  Although  the  various  contributors  have  not  given  an  exhausti.c 
account  of  their  subjects,  it  is  believed  that  there  is  sufficient  information  for  the  reader  to 
evaluate  the  tradeoffs  involved,  for  example,  in  choosing  a  particular  mode  of  marking  a 
given  location  or  devising  an  underwater  navigation  system. 

Section  !  is  a  series  of  charts  which  present  the  tradeoffs  involved  in  various  methods 
of  underwater  information  transfer  or  site  marking.  (The  charts  are  ir.  an  extremely  con¬ 
densed  form.)  These  charts  are  self-contained,  and  can  be  used  without  referring  to  their 
corresponding  texts  in  Section  II.  It  is  strongly  recommended,  however,  that  the  appropriate 
text  be  consulted  because  some  assumptions  made  by  the  authors  might  vary  from  the 
specific  parameters  envisioned  for  the  reader's  sy  stem.  The  charts  particularly  emphasize 
maximum  range  for  threshold  detectability  because  the  background  noise,  dispersion, 
and  attenuation  encountered  in  the  ocean  often  place  severe  restrictions  on  the  usable 
range.  The  ranges  on  the  charts  reflect  as  much  as  possible  the  natural  environmental 
limitations  on  range.  Considerations  of  additional  noise  created  by  detection  or  h\  inMiffi- 
cient  state-of-the-art  detection  sensitivity  are  discussed  in  Section  11. 

Section  II.  which  contains  articles  corresponding  to  the  charts  in  Section  1.  explains 
in  detail  the  methods  of  calculating  the  values  which  appear  on  the  charts.  In  addition,  the 
problems  involved  in  implementing  the  various  methods  are  iiscussed.  In  view  of  the 
sponsor's  interest  in  navigation  applications,  some  of  the  articles  are  written  more  directly 
towards  that  goal,  although  much  of  the  material  is  applicable  in  other  areas  The  last 
article  in  Section  II.  necessarily  more  qualitative  than  the  others  and  lacking  a  correspond¬ 
ing  chart,  reviews  undersea  mechanical  navigation  aids. 


SECTION  I:  CHARTS 


This  section  consists  of  charts  which  tabuiate  the  tradeoffs  between  the  various 
methods  of  underwater  marking  and  information  transfer.  The  juxtaposition  of  the  charts 
makes  comparisons  convenient.  The  charts,  however,  should  be  used  only  in  conjunction 
With  their  corresponding  articles  in  Section  II  because  the  charts  are  the  direct  contributions 
of  the  various  authors  and  are  written  from  slightly  different  perspectives. 

Some  charts  were  prepared  assuming  coherent  detection,  while  others  assumed 
incoherent  detection.  Slight  differences  in  the  signal-to-noise-retio  criteria  appear  among 
the  various  authors.  The  ranges  in  some  of  the  charts  are  obtainable  with  current  -  quipment 
or  techniques,  while  others  are  based  on  what  is  available  in  the  near  term,  ror  chest  re  a 
sons,  the  charts  must  be  used  with  care,  and  the  corresponding  articles  should  oe  read  before 
.tcvsiingful  comparisons  are  made. 
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STATIC  ELECTRIC  AND  MAGNETIC  FIELDS 


INTRODUCTION 

A  technique  that  has  been  used  to  detect  and  locate  submerged  objects  is  the  observa¬ 
tion  of  their  associated  electric  or  magnetic  fields.  Objects  consisting  of  more  than  one  type 
of  metal  usually  experience  electrolytic  decomposition  in  the  sea.  The  electric  field  caused 
by  this  process  is  detectable  at  a  distance  generally  greater  than  the  dimensions  of  the  object 
being  sought.  These  fields  and  the  objects  producing  them  have  been  located  by  tewing 
electrodes  from  surface  vessels,  a  method  which  has  been  used  to  locate  sunken  torpedoes, 
ships,  and  planes  (ref.  1 ). 

Towed  magnetometers  nave  also  been  used  to  locate  submerged  objects  which  are 
magnetized  or  magnetizable  in  the  earth’s  field.  This  type  of  object  superimposes  on  the 
earth's  field  a  short  wavelength  anomaly  that  can  be  detected  by  the  magnetometer.  The 
magnetometer  trace  of  the  sunken  USS  THRESHER  is  a  good  example  (ref.  2). 

The  present  article  is  concerned  with  the  detection  of  various  types  of  static  electric 
and  magnetic  sources  which  are  located  on  or  near  the  ocean  bottom.  It  will  be  assumed  that 
these  objects  arc  bei”?  searched  for  with  a  detector  that  is  also  near  the  bottom,  for  example, 
a  magnetometer  mounted  on  a  towed  fish  or  submersible.  The  basic  questions  that  this 
article  will  address  are  listed  below. 

1.  What  is  the  maximum  range  at  which  an  electric  or  magnetic  source  can  be 
detected  given  a  set  of  parameters  describing  the  source  and  surrounding  environment? 

2.  How-  does  the  maximum  range  change  as  these  parameters  arc  varied? 

In  the  calculations  which  follow,  two  assumptions  ensure  that  the  ranges  obtained 
are  the  upper  bounds  to  ranges  which  can  be  obtained  in  practice.  In  the  first  condition, 
both  the  signal  and  the  noise  statistics  are  assumed  known  and  coherent  detection  is  used. 

In  an  actual  attempt  to  locate  an  object,  often  very  little  can  be  assumed  about  the  expected 
type  of  signal  and  noise.  If  bounds  can  be  placed  on  the  distance  and-or  the  direction  to  the 
souree.  however,  some  information,  which  can  be  used  in  processing  the  received  signals,  is 
available  concerning  the  signal's  amplitude,  frequency  spectrum,  etc.  The  greater  the 
a  priori  knowledge  of  the  expected  signal,  the  closer  the  ranges  obtained  in  an  actual  search 
will  approximate  those  given  in  this  article.  The  second  factor,  which  probably  establishes 
the  ranges  as  the  upper  bounds,  is  the  neglect  of  any  increase  in  ambient  noise  due  to  the 
detection  process.  Possible  noise,  caused  by  equipment  aboard  the  submersible,  turbulence 
created  by  the  platform,  electrolytic  action.  etc_  h3S  not  been  considered.  Only  the  noise 
processes  over  which  the  searcher  has  no  control  have  been  included  in  the  present 
calculations. 

Certain  parameters  pertaining  to  the  ocean  environment  and  to  the  method  of  search 
have  been  assigned  numerical  values  which  will  be  listed  in  this  article.  The  dependence  of 
the  detection  range  on  some  of  these  parameters  is  given  in  Chart  1 .  Section  I.  Both  the 
source  and  detector  are  assumed  near  the  bottom  ir.  4  km  of  water;  this  depth  is  almost  the 
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average  depth  of  the  oceans  (ref.  3).  The  static  magnetic  field,  Bq,  has  been  assigned  a  value 
of  40.000  gammas,  and  the  seawater  conductivity.  <7,  is  4  mhos/m.  The  bottom  current  and 
the  velocity  of  the  detector  are  0.1  knot  and  1.0  knot,  respectively.  The  permeability  of  the 
ocean  and  the  seabed  have  been  assigned  the  value  of  the  permeability  of  free  space. 

|x=  =  4«  x  10-/  N/A-.  The  rationalized  mks  system  is  used  throughout  this  article. 

BASIC  APPROACH 

The  search  scheme  used  is  shown  in  figure  1 . 


s 


leurc  I.  Search  kbraie. 


Tlie  platform.  P.  carrying  the  detector  is  assumed  to  move  along  a  straight  horizontal 
line  with  the  detector  at  the  same  depth  as  the  source.  S.  The  platform  is  assumed  to  start  at 
some  point. -Xq- along  this  line.  It  then  proceeds  with  the  detector  past  the  point  of  closest 
approach  of  the  search  path  with  the  source  and  on  to  X  =  +Xq-  At  this  point,  a  decision  is 
made  concerning  the  presence  or  absence  of  the  source.  The  orientation  of  the  source's 
dipole  moment  is  assumed  horizontal  and  perpendicular  to  the  platform's  path.  The  angle 
between  the  dipole  moment  of  the  source  and  the  line  drawn  from  the  source  to  the  detector 
is  0.  The  distance  along  the  platform's  path  measured  from  the  point  of  closest  approach. 

X  =  0.  is  then  given  bv  X  =  R  tan  0.  In  the  following  discussion,  we  will  be  concerned  with  the 
spatial  variation  of  signals  and  the  sparial  spectra!  Jensilv  of  noise,  although  one  can  also 
consider  their  time  variation,  as  seen  bv  the  detector,  with  equivalent  results.  The  connection 
between  the  lime  and  the  spatial  variation  is  made  through  the  relation.  X  =  YT  Xq.  where 
V  is  the  platform's  velocity  and  T  is  the  time. 

The  detector  is  assumed  to  be  a  single  component  (scalar)  detector,  such  as  a  flux- 
gate  magnetometer  or  electric  dipole.  It  is  assumed  that  this  detector  is  alwavs  oriented  so 
that  tne  signal  lies  along  its  sensitive  direction.  This  is  consistent  w  ith  our  assumption  of  a 
coherent  search.  We  will  also  assume  that  the  noise  is  independent  of  the  antenna's 
orientation. 

The  field  due  to  an  infinitesimal  dipole  in  the  infinite  ocean  is 

F  = - D  —  12cos0f*sin  0#).  (1) 

(x-  R-r  ~ 

where  D  =  p  4ztrwhen  F^  =  IT.  the  electric  field,  and  D  =  iiqM/4it  when  F  =  B.  the  magnetic 
field  (refs.  4  and  5).  Tne  current  and  magnetic  dipole  moments  of  the  sources  are  p  and  M. 


respectively.  The  unit  vector,  r.  is  directed  from  the  source  to  the  detector,  and  0  is  a  unit 
vector  pointing  in  the  direction  of  increasing  0.  The  use  of  Eq.  1  is  justified  for  the  sources 
considered,  providing  we  are  far  from  the  source  compared  with  its  dimensions.  Because  both 
the  source  and  detector  are  near  the  bottom,  the  electric  field  for  a  constant  dipole  moment 
can  rise  to  a  value  twice  that  given  by  Eq.  3  (ref.  5).  In  view  of  the  porosity  of  the  seabed 
and  the  inherent  uncertainties  in  the  ambient  noise.  Eq.  1  will  be  used.  For  the  magnetic 
dipole.  Eq.  1  is  correct  provided  that  the  permeability  of  the  ocean  and  the  seabed  are  equal 
(as  assumed)  and  that  the  coil  is  sufficiently  insulated  so  that  its  electric  field  is  negligible 
outside  the  coil. 

At  the  end  of  the  tow.  at  point  Xq.u  decision  is  made  concerning  the  presence  or 
absence  of  the  source.  This  decision  is  effected  by  passing  the  output  of  the  detector 
through  the  appropriate  linear  filter,  which  is  matched  to  the  signal  and  ambient  noise. 

This  “matched  filter"  technique  will  maximize  the  filter's  output  signal-to-noise  ratio  at  point 
XQ(rcfs.  6  and  7).  The  output  at  point  X^is  then  compared  with  a  fixed  decision  level, 
related  to  the  detection  and  false  alarm  probabilities,  and  the  presence  or  absence  ot  the 
source's  signal  is  decided. 

For  an  infinite  measuring  time,  the  relationship  between  the  signal,  noise,  and 
signal-to-noise  ratio  ( ref.  S )  is  given  by 


d?  = 


iSi(u)j- 

4>(cj) 


<2> 


where  d  is  the  signal-to-noise  ratio,  u  is  the  time  radian  frequency .  Stcu)  is  the  time  Founer 
transform  of  the  signal,  and^tcj)  is  the  lime  noise  spectral  power  densitv .  Using  the  relation¬ 
ships  «  =  VK.  |SM-  =  1|S(K)|-I  ’IJVj-l.  and<p(K>=  Eq.  2  becomes 


where  K  is  the  spatial  angular  frequency.  V  is  the  detector’s  velocitv .  St  K  >  is  the  spatial 
Fourier  transform  of  the  magnitude  of  Eq.  I.  and  3>(K)  is  the  spatial  noise  spectral  power 
density.  It  is  interesting  to  note  that  the  form  of  Eqs  2  and  3  is  indifferent  as  to  whether  we 
measure  in  the  time  or  spatial  domain:  however,  wx-  will  continue  to  use  the  spatial  viewpoint 
because  wx*  arc  dealing  with  static  sources.  Two  assumptions  implicit  in  Eqs.  2  and  3  should 
be  mentioned.  Both  Eqs.  2  and  3  refer  to  infinite  tow  paths:  however,  because  of  the  rapid 
decay  of  Eq.  I  for  iaige  X.  Eqs.  2  and  3  apply  providing  wx*  take  j.Xpj  >3R.  Ai>o.  both 
Eqs.  2  and  3  are  applicable  only  to  noise  which  is  stationary  and  gaussian. 

Taking  the  Fourier  transform  of  Eq.  1 .  we  find  that  jS(  K  )j-  is  almost  constant  for 
jK|  <  1  R  and  decreases  rapidly  for  }K|  >  l  /R.  In  view  of  this  property,  wx*  take  IS(K)j-  to  5>e 
constant  for  |K|  <  1  'R  and  zero  otherwise.  With  this  simplification  and  the  use  of  ParsevaEs 
theorem,  we  obtain 


/*00 

*  df 

f  n?!^x 

— oO 

0 

$(K) 

f«  J 

(4) 


where  fu  is  the  effective  upper  cut-off  frequency: 

f  =ik-_L 

U  In  2-R  ' 

The  left-hand  integral  in  Eq.  4  has  the  value 

OO 

I  [f!2  dx  =  D23  (2v)6 

-OO  * 

where  r  =  2»  R  and  the  term  in  the  bracket  in  Eq.  4  is  the  average  of  the  reciprocal  of  the 
noise  density  over  the  spectrum  of  the  signal,  i.e..  0  <  f  <  l/2oR. 

Both  the  magnetic  and  electric  noise  were  fit  to  a  noise  power  density  of  the  form 

<f»e.  m  (K)  =  — e-rm-  +  Ce.  m.  (5B) 

k 

where  e  refers  to  the  electric  noise,  m  refers  to  the  magnetic  noise,  and  Be.m  and  Ce.m  are  a 
set  of  four  constants.  For  this  noise  spectrum,  the  brackcicd  term  in  Eq.  4  is 


Eqs.  5  and  6  can  oe  combined  to  give 


where  G  =  D2/C  (2s-)6  (21/32). 

Hie  value  chosen  for  d.  the  signal-to-noise  ratio,  will  determine  the  false  alarm  proba¬ 
bility  and  the  detection  probability  for  a  given  run.  In  the  present  calculations,  a  value  of  two 
is  used.  This  value  corresponds  to  a  detection  probability  of  0.9  and  a  false  alarm  probability 
of  0.2  (ref  S).  For  the  same  value  of  d.  one  can  obtain  a  greater  probability  of  detection  at 
the  expense  of  a  greater  false  alarm  probability  and  vice  versa.  Greater  values  of  d  will  enhance 
the  probability  of  detection  or  lessen  the  false  alarm  probability  when  the  other  is  constant. 
The  dependence  of  the  detection  range  on  d  is  given  in  Chart  1 .  Section  I. 

Before  using  Eq.  7  to  calculate  numerical  results  for  the  electric  and  magnetic  cases, 
we  will  derive  its  asymptotic  forms.  For  r«(2sC)/B.  (Br)/(2«C)  «  1  and 

/  C  \'5 

^[pj  ■  ,s> 

For  r  » (2rC)/B.  (rB)/(2«C)  »  1  and 
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1/6 


(9) 


Thus,  r  varies  as  the  one-fifth  power  of  G  for  small  r  and  as  the  one-sixth  power  of  G 
for  large  r,  where  r  is  to  be  compared  with  tq  =  (2~C)/B. 

MAGNETIC  FIELDS 


The  first  type  of  magnetic  source  to  be  considered  is  a  copper  torus  through  which 
current  is  passed.  The  coil  is  assumed  sufficiently  insulated  to  ensure  that  there  is  essentially 
no  current  outside  the  windings.  The  magnetic  moment.  M.  of  such  a  coil  is  given  by 


where  A  is  the  cross-sectional  area  (through  the  windings).  P  is  the  power  input,  r  j  is  the 
large  radius  of  the  torus,  and  p  is  the  resistivity  of  the  windings.  Using  values  of  A  =  "(0.1  m)~ 
=  3.14  x  10"-  m2.  P=  10W,  rj  =  1  m,  and  p  =  1.7  x  I0~8o-m  (j.e..  the  resistivity  of  copper), 
one  obtains  a  magnetic  moment  of  5.44  x  1(P  A-m^.  The  magnetic  moment  for  different 
types  of  coils  and  input  power  can  be  calculated  from  Eq.  10.  It  is  interesting  to  note  that 
larger  magnetic  moments  for  a  fixed  coil  displacement  arc  possible  by  increasing  the  coil's 
diameter  at  file  expense  of  its  cross-sectional  area.  For  the  parameters  chosen,  the  coil's 
weight  is  about  2  tons  and  displaces  0.2  m-7.  For  a  30-day  operating  life,  this  coil  will 
require  170  lb  of  silver-zinc  batteries  which  displace  about  2  ftJ. 

The  second  type  of  source  to  be  considered  is  the  permanent  magnet.  This  source 
has  the  advantage  of  not  requiring  any  power,  but  the  disadvantage  of  not  haring  its  magne¬ 
tization  readily  controllable  once  it  is  in  the  w-ater.  Using  Alnico  5.  we  can  obtain  a  perma¬ 
nent  magnet  with  a  remnant  magnetic  induction  of  13.100  G  (ref.  9).  which  corresponds  to 
a  magnetization  of  1.05  x  JO*7  A/m.  Again  assuming  a  magnet  weight  of  2  tons,  we  obtain 
a  volume  of  0.25  nv7.  yielding  a  net  dipole  moment  of  2.61  x  lO^A-m-  for  this  source. 

A  third  type  of  source  is  a  soienoid  th3l  is  filled  with  ?  high  permeability  material. 

Tire  magnetic-field  intensity  in  a  solenoid,  which  is  long  compared  with  its  diameter,  can  be 
written  as  H  =(  1  /(2rr^L)l  (PW/p  ^  A,)I/2„  where  r-  is  the  solenoid's  radius.  L  is  its  length. 

W  is  the  weight  of  the  windings,  and  p  and  pw  are  the  winding's  resistivity  and  weight  density, 
respectively.  This  formula  is  strictly  true  only  as  the  length-to-diame:er  ratio  approaches 
infinity;  however,  because  of  the  large  value  of  the  ratio  in  the  present  case,  it  is  an  excellent 
approximation.  Choosing  values  of  =  0. 1  m.  L  =  2~  m.  P  =  70W.  and  W  =  2  tons,  the  result¬ 
ing  field  intensity  for  copper  (pw  =  S.96  x  lCPkg/m^.  p  =  1 .7  x  1  CT^n-m)  is  ~93  Oe.  Filling 
the  solenoid  with  a  high  permeability  material  will  lower  the  m3gnetic-field  intensity  from 
the  value  previously  given  due  to  the  opposing  magnetization  field.  For  this  geometry,  the 
demagnetization  factor.  N/4*.  is  4.3 1  x  1  O'-7  (ref.  9).  .Assuming  the  solenoid  to  be  filled 
with  Supemrendur.  a  net  magnet ic-fieJd  intensity  of  10  Oe  is  obtained  from  its  magnetization 
curve  (ref.  9)  and  the  demagnetization  factor.  At  this  field  strength,  the  magnetic  induction 
is  20.000  G.  and  the  magnetization  is  1.6  x  I0*7  A/m.  which  is  near  saturation.  The  perme¬ 
able  core  has  a  volume  of  ~0.2  mJ  so  that  the  total  magnetic  dipole  moment  is  3.2  x  1 0^  A-m 
In  addition  to  the  2-ton  weight  of  the  copper,  the  Supemrendur  weighs  about  l.S  tons  and  the 
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batten-  weight  for  1  month  of  operation  is  about  0.6  toil.  It  is  interesting  to  note  that  the 
copper  coil  previously  considered,  w'hcn  run  at  a  power  level  of  70  W,  will  yield  a  dipole 
moment  of  1 .4  x  10^  A/m,  a  factor  approximately  20  times  lower  than  that  for  the  filled 
solenoid. 

The  last  source  to  be  considered  is  a  superconducting  coil.  Recent  proposals  have 
indicated  the  feasibility  of  constructing  large  superconducting  coils  with  magnetic  moments 
of  approximately  10'  A-m^.  For  a  superconducting  coil  of  1-m  radius,  a  magnetic  moment 
of  3.56x  106  A  -m-  appears  feasible  using  a  matrix  of  type  II  superconductor  and  normal 
metal.  Because  the  critical  field  for  the  best  superconductors  is  about  one  order  of  magnitude 
higher  than  that  for  the  saturation  induction  of  the  highest  permeability  materials,  the  order- 
of-magnitude  increase  in  the  dipole  moment  upon  moving  from  the  filled  solenoid  to  the 
superconducting  coil  is  reasonable. 

Natural  magnetic  noise  arises  from  such  temporal  sources  as  solar  activity,  ionospheric 
activity,  and  the  motion  of  seawater  through  e  earth's  magnetic  field.  Except  for  occa¬ 
sional  disturbed  periods  and  the  diurnal  variations,  these  noise  sources  are  associated  with  dailv 
excursions  of  approximately  a  few  gammas  { 17  =  1  O'3  G) .  In  a  search  for  static  magnetic 
fields.  Spiess.  Mudie.  and  Lowenstien  (SML)  (ref.  10)  have  shown  that  the  major  environ¬ 
mental  limitation  comes  from  short  wavelength  anomalies  superimposed  on  the  earth's  dipole 
field.  Their  results  will  be  used  to  estimate  the  magnetic  noise  that  might  be  encountered  on 
a  search  as  proposed  in  this  article.  Due  to  the  necessity  of  interpolating  data  from  SML's 
graphs,  an  error  in  our  noise  figures  of  a  factor  of  two  should  be  allowed.  Because  the  range 
varies  as  the  fifth  power  of  the  reciprocal  of  the  average  noise,  this  error,  however,  represents 
only  a  small  error  in  the  estimated  range. 

Based  on  Eq.  5B.  two  relationships  will  be  needed  to  specify  the  noise  power  density. 
Thus,  an  expanded  scale  version  of  SML’s  magnetometer  trace  was  obtained  and  a  10-km 
portion  was  digitized  and  Fourier  analyzed.  It  was  noted  that  the  resulting  spectrum  increased 
rapidly  for  spatial  frequencies  less  than  0.001  cpm.  Because  in  Eq.  5B  increases  rapidlv  for 
f  <  B/(  2ttC).  B.'(2«C)  was  set  equal  to  0.001  Hz.  A  second  relationship  was  obtained  by  relat¬ 
ing  the  varianc*  of  the  filtered  spectrum  given  by  SML  to  the  integral  of  the  noise  power  spec¬ 
trum  across  their  bandpass.  The  rms  value  of  the  filtered  noise,  obtained  from  ihe  SML 
spectrum,  was  halved  in  an  attempt  to  include  the  fact  that  their  spectrum  was  taken  in  an 
area  of  "high  magnetic  roughness. "  The  resulting  rms  noise  was  thm  reduced  bv  another 
factor  of  I.V3  to  account  for  the  reduction  in  noise,  assuming  isotropy,  that  would  be 
obtained,  by  replacing  the  SML  total  field  magnetometer  with  the  one-component  magne¬ 
tometer  assumed  in  this  article.  Tnis  final  adjusted  value  of  nns  noise  was  compared  with  the 
integral  of  4> across  the  SML  bandpass  to  \ Icid  the  second  condition  on  B  and  C.  The  result- 
ing  values  arc  B  =  49.627-  and  C  =  7.897  x  1CP  Y-m.  These  values  plus  the  values  for  the 
magnetization  and  the  formulas  previously  given  are  sufficient  to  estimate  the  maximum  range 
in  a  magnetic  search. 

ELECTRIC  FIELDS 

For  the  present  calculation,  the  source  of  the  electric  field  is  considered  to  be  an 
insulated  linear  conductor  terminated  in  conducting  spheres.  It  is  assumed  that  the  diameter 
of  the  conductor-insulator  combination  is  small  when  compared  with  the  diameter  of  the 
sphere  and  that  the  IR  drop  along  the  wire  is  negligible  when  compared  with  the  applied 
FMF.  The  EMF  is  assumed  to  be  inserted  along  a  length  of  the  antenna  and  to  have  an  inter¬ 
nal  resistance  that  is  small  compared  with  the  antenna's. 
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Because  of  the  remarks  made  in  the  section,  Basic  Approach,  the  input  resistance  will 
be  that  of  the  antenna  in  infinite  space.  Because  of  the  assumption  made  concerning  the 
antenna  wire,  the  input  resistance  will  be  the  resistance  between  the  two  terminating  spheres 
in  the  absence  of  wire.  For  this  case,  we  can  use  the  resistance  capacitance  analogs'  which 
relates  the  resistance,  R,  between  the  spheres  in  an  infinite  medium  of  conductivity,  cr.  to 
their  capacitance.  C.  in  vacuo  through  the  relationship  R  =  tQ/crC’fref.  1 1 ).  The  capacity 
between  two  spheres  of  radius  r$  with  their  centers  a  distance  d  apart  is  given  by 

OG 

O  =  2~CQrs sinh P  Icsch(2n  -  1  )&  +  csch  2n/Jj . 
n=l 

where  cosh/J  =  d/2r$  (ref.  1 2).  Choosing  d  =  2«  m  and  rs  =  0.1  m,  we  have  d  »rs  for  which 
C*  =  2-cprs.  The  resistance  is  then  given  by 
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From  Eqs.  5  and  7. 


4(444). 

C  y*trc t]  \>2/ 


where  p  =  Id  =  (P/R)^'~  (d  is  the  dipole  moment).  Using  values  of  r$  =  0.1  m  and  d  =  2»  m. 
we  find  that 

1.5S5  x  103  |p|V--m4 


where  |P|  is  the  power  expressed  in  watts. 

Natural  electric  noise  arises  from  three  sources.  The  first  source  corresponds  to  tem¬ 
poral  fluctuations  arising  outside  of  the  earth  and  propagating  beneath  the  Oe can’s  surface. 
Using  noise  figures  given  by  Liebermann  (ref.  1 3)  and  Cox  ei  al  (ref.  14)  and  assuming  an 
infinite!)  deep  sea.  we  find  the  electric  noise  ranges  from  1.0  to  4.0  \  10'*-*  fV  m)-'Hz  for 
frequencies  between  0.003  to  0.001  Hz.  Viewed  as  a  spatial  variation,  this  corresponds  to  a 
noise  power  density  of  0.51  to  2.6  x  10"'  4  (V/m)-/cpm  for  frequencies  between  0.002  and 
0.006  cpm  for  a  detector  speed  of  I  knot. 

The  interaction  between  the  movement  of  the  water  and  the  ambient  magnetic  field 
will  also  cause  an  electric  field.  The  magnitude  of  this  field  depends  to  some  extent  upon 
the  ability  of  the  surrounding  water  to  "short  circuit"  the  potential  across  the  water  in  the 
vicinity  of  the  receiving  antenna.  For  the  complete!)  short-circuited  case,  the  electric  field 
along  the  antenna  is  given  by  the  product  of  Vjj  and  By.  where  Vjj  is  the  component  of 
velocitv  perpendicular  to  the  antenna  and  By  is  the  vertical  magnetic  field  (ref.  15).  If  we 
divide  both  By  and  V^j  into  their  (spatiallv )  constant  and  varving  components,  ‘.he  resulting 
product  will  contain  four  terms.  One  of  these  terms  will  be  constant  3nd  can  be  neglected 
because  it  will  be  rejected  by  our  matched  filter.  A  second  term,  relating  the  two  spatiallv 
van  ing  parts,  will  be  small.  The  remaining  two  terms  represent  the  interaction  of  the  constant 
part  of  the  velocity  with  the  varying  part  of  the  magnetic  field  and  vice  versa.  Tire  former 
term  is  the  product  of  the  square  of  the  mean  velocitv  and  the  magnetic-noise  power  spectrum. 
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which  was  calculated  earlier.  At  a  frequency  of  0.001  cpm,  this  is  4  x  10"*  5  (V/m)-/cpm. 

Both  these  sources  are  dominated  by  velocity  fluctuations  interacting  with  the  earth’s 
field.  The  velocity  fluctuations  used  in  this  article  were  measured  0.75  to  1.75  m  off  the  bot¬ 
tom  in  water  1 2  to  20  m  deep  in  an  area  off  Red  Wharf  Bay  (~53°20’N  -  4°  1 0'W).  Whether 
the  form  of  the  turbulence  spectrum  is  maintained  in  deeper  water  has  not  been  determined. 
Four  points  of  the  low-frequency  end  of  the  Red  Wharf  Bay  data  were  fit  to  a  straight  line  of 
the  form  $(cu)  =  A’/fw  +  B\  where  $(cj)  is  the  turbulent  noise  power  measured  on  the  bottom. 
fw  is  the  time  frequency,  and  A’  and  B’  are  constants  depending  on  the  mean  current.  To 
obtain  the  needed  values,  it  was  necessary  to  extrapolate  the  observed  data  by  a  factor  of  1 0: 
however,  the  straight  line  fit  had  a  correlation  coefficient  of  0.97.  The  velocity  noise  spec- 
trum.  viewed  as  carried  along  with  the  mean  current,  was  then  found  to  be$(K)  =  (0.102U)- 
(7.93Usec  +  0.04656/f$).  where  U  is  the  rms  current  and  fs  is  the  spatial  frequency.  The  elec¬ 
tric  noise  power  density  will  be  Bv~  times  this  value.  Inserting  values  of  U  =  0.1  knot  and 
8V  =  40.000  gammas  and  comparing  the  result  with  Eq.  5B,  we  obtain  B  =  1 .3  x  1  O'*4  (V/m)~ 
and  C  =  1.8  x  10~*4  m  ( V/m )-  for  the  electric-noise  coefficients.  These  coefficients,  plus  the 
expression  for  the  electric  signal,  allow  us  to  calculate  range  as  a  function  of  input  power  by 
the  method  outlined  in  the  section,  Basic  Approach. 

It  should  be  noted  that  electric  “velocity"  noise  dominates  “atmospheric"  noise  only 
if  the  mean  current  exceeds  a  certain  minimum.  For  the  values  of  noise  in  this  article,  this 
minimum  current  is  0.1  cm/sec.  The  present  analysis  does  not  extend  to  currents  lower  than 
this  value. 

RESULTS  aND  DISCUSSION 

Figure  2  shows  the  variation  of  range  with  increasing  magnetic  moment.  It  is  interest¬ 
ing  to  note  that  increasing  the  source's  magnetic  moment  by  a  factor  of  HP  increases  the 
maximum  range  by  less  than  a  factor  of  10.  Near  a  magnetization  M  =  1  x  10-  A-nr*.  the 
range  increases  as  R  ~  mB-39G  •vhiic  for  M  near  1x10^  A-m-.  R  ~  .\f0  348  Over  the  entire 
range,  i  x  10-  A-m-<  M  <1x10^  A-m-.  R~  The  variation  of  the  slope  with  mag¬ 

netization  is  indicative  of  the  competition  between  the  “white”  and  1/f  portions  of  the  noise 
spectrum.  Knowing  the  variation  of  R  with  M  and  the  dependence  of  M  on  various  parameters, 
the  variation  of  R  with  these  parameters  can  be  determined.  This  dependence  is  listed  in 
Chart  1.  Section  I.  The  variation  of  R  with  these  parameters  is  unexceptional  except  for  the 
large  increase  in  range  with  the  diameter  of  the  copper  coil  as  compared  with  an  increase  in 
its  input  power. 

The  range  of  an  electric  dipole  3S  a  function  of  input  power  is  given  in  figure  3.  The 
range  obeys  dependence  on  the  power  at  all  power  levels  between  10'*  and  1 04  W. 

This  behavior  is  due  to  the  dominance  of  I/f-type  noise. 

It  is  interesting  to  note  the  strength  of  the  signal  at  the  distance  of  closest  approach 
for  the  magnetic  and  electric  cases.  At  this  point,  the  magnetic  field  has  a  value  of  74  gammas 
fora  dipole  strength  of  100  A-m-  and  a  value  of  22  gammas  for  a  dipole  strength  of  107  A-m-. 
These  values  arc  well  within  th-  ocnsitivity  of  conventional  magnetometers,  and  sensitivities  of 
this  type  have  been  exhibited  in  the  ocean  environment.  For  the  electric  fieid.  the  potential 
across  the  ends  of  an  antenna  of  length  2«  m  will  be  about  2  p.V.  This  level  can  be  measured 
by  standard  potentiometer  techniques:  however,  care  must  be  exercised  to  minimize  thermal 
EMF’s.  turbulence  introduced  into  the  water  by  the  detector  or  platform,  and  electrolytic 
action  on  the  receiving  electrodes. 
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The  question  still  remains  as  to  whether  the  ranges  predicted  in  this  article  are  typical 
of  the  oceanic  environment  as  a  whole:  Range  will  change  in  different  geographical  locations 
due  to  variations  in  the  effective  noise  power  density,  if  it  is  assumed  that  the  form  of  this 
function  changes  very  little  from  the  form  in  this  article,  the  noise  spectrum  will  be  proportional 
to  the  variance  of  the  noise  across  the  passband.  For  the  magnetic  case,  it  is  estimated  that  on 
a  worldwide  basis  the  rms  noise  can  vary  by  a  factor  o‘  six;  this  implies  that  the  obtainable 
range  can  vary  by  a  factor  of  about  two.  For  the  electrical  case,  deep  ocean  currents  can 
assume  values  of  approximately  O.Oi  to  0.5  knot  (ref.  16).  These  values  imply  ranges  varying 
from  a  factor  2.3  times  greater  to  a  factor  J  .8  limes  smaller  than  those  in  this  article.  The 
large  variability  of  the  maximum  range  obtained  thus  emphasizes  the  dependence  of  search 
effectiveness  on  environmental  parameters. 
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DETECTION  OF  RADIOACTIVE  BEACONS  IN  THE  OCEAN 


i 


Two  approaches  to  the  detection  of  radioactive  beacons  in  water  are  discussed  in  this 
article:  detection  of  a  point  source  of  radioactivity  on  the  bottom  or  in  open  water  and  detec¬ 
tion  of  induced  activity  from  a  neutron  source  in  open  water. 

DETECTION  Or  A  RADIOACTIVE  POINT  SOURCE  UNDERWATER 


The  range  of  nuclear  radiation  in  water  is  quite  limited.  This  is  true  whether  consider¬ 
ing  alpha,  beta,  or  gamma  rays,  or  neutrons,  all  of  which  are  readily  attenuated  b>  water. 

For  gamma  rays,  the  ntensity  from  a  point  source  can  be  calculated  from  the  following 
relationship: 


1  = 


Nr 


Bc~FR 

4rR~ 


where  R  is  the  range  from  the  source:  pis  the  linear  attenuation  coefficient  (for  energy):  B  is  a 
buildup  factor  to  account  for  photons  that  have  been  scattered  and  lost  energy,  but  have  not 
been  lost  to  the  beam:  Nq  is  the  original  intensity:  and  1  is  the  photon  density  at  R.  At  very 
short  ranges  (pR  «  1 ).  the  inverse  square  term.  l/(4rrR-).  controls  the  change  in  photon  den¬ 
sity  as  a  function  of  distance.  At  three  inean-free-paths  (the  distance  Rq  =  Ip),  the  inverse 
square  and  attenuation  are  equally  effective:  beyond  three  Rq.  the  attenuation  term  becomes 
progressively  more  dominant.  In  air.  Rq  is  approximately  hundreds  of  feet  (419  ft  at  *.0  MeV) 
(ref.  1 ).  but  underwater  Rq  is  very  short  ( 14  cm  at  1  McY  and  a  maximum  of  60  cm  for  energies 
greater  than  25  MeV).  For  a  100-Ci  point  source  of  6^Co.  the  detection  range  underwater  is 
approximately  9  ft  (1  is  equal  to  the  gamma-ray  background  of  seawater).  Even  with  a  ingh- 
energy  machine  (25  MeV).  the  range  for  a  yield  equivalent  to  100  Ci  will  be  no  more  than  2S  ft. 
At  best,  gamma-ray  sources  are  short-range  beacons. 

Neutron  sources  are  equally  short  range  when  considering  detection  of  the  source  neu¬ 
tron.  For  fast  neutrons  from  a  fission  source,  i.e..  -^-Cf.  the  flux  density  decreases  by  3  .actor 
of  lOevery  10cm  (3l  a  1-ft  range)  to  20cm  (at  a  3-ft  increase  in  range).  Neutrons  from  a  10*  * 
n/sec  source  (a  large  source  by  present  standards)  will  hardly  be  detectable  bey  ond  8  ft. 

Among  the  nuclear  particles,  only  neutrinos  penetrate  water  readily .  Their  low  inter¬ 
action  rate  with  matter  makes  detection  extremely  inefficient  and  establishes  a  very  high 
background  flux,  primarily  from  extraterrestial  sources.  Both  facets,  the  lack  of  a  suitable 
detector  and  the  probable  high  background,  negate  consideration  of  a  neutrino  beacon. 

DETECTION  OF  NEUTRON-INDUCED  ACTIVITY  IN  SEAWATER 


Most  neutrons  from  a  point  source  in  water  are  themialized  and  captured  within 
approximately  1  ft  of  the  source.  In  the  process.  3  scries  of  radioactive  isotopes  is  induced 
in  the  seawater  surrounding  the  source.  Most  prominent  among  the  gamma-ray  emitters 


induced  by  this  process  are  J^C1  (38  min)  and  -^Na  (15  hr).  If  the  seawater  surrounding  the 
source  flows  past  the  source,  the  induced  activity  wfll  be  carried  from  the  source  in  the  same 
manner  that  smoke  is  carried  by  air  from  a  campfire.  The  induced  activity  can  be  detected,  at 
a  considerable  range  from  a  point  source,  as  a  distributed  source  of  active  water.  For  a  neutron 
source  of  10' '  n/sec  in  a  current  of  1  knot,  the  induced  activity  in  seawater  passing  within  1  ft 
of  the  source,  is  approximately  6x1 0"^  Ci/ 1 .  This  is  2  x  10-  times  the  natural  background  of 
seawater.  Assume  that  the  cross  sectional  area  of  this  stream  of  active  solution,  which  is  approxi¬ 
mately  0.28  m-.  gradually  increases  as  it  «s  carried  away  by  the  current.  The  activity.  1  hr 
later,  will  be  1  km  downstream  with  a  frontal  area  of  500  m-.  a  nominal  diffusion  growth 
(ref.  2):  however,  it  will  still  be  readily  and  unambiguously  detected  with  present  underwater 
detection  systems.  If  the  currents  are  negligible  and  if  the  activity  drifts  around  so  that  the 
general  region  becomes  activated  (i.e..  smokey).  10*  *  n/sec  will  induce  enough  activity  to 
double  the  background  in  8  x  10-*  m-*  of  seawater.  In  practical  situations,  the  activity  will 
not  be  uniformly  distributed  and  the  “smoke”  will  constitute  a  trail  to  be  followed  to  the 
source.  In  the  final  approach  to  a  neutron  source  in  seawater,  the  -aptured  gamma  rays, 
which  are  emitted  each  time  a  neutron  is  captured,  will  clearly  define  the  area  (ref.  3).  These 
captured  rays  from  seawater  are  high-energy  gamma  rays  that  can  probably  be  detected  at 
16-f!  '-'>ges  from  a  10'  '-n/sec  source. 

The  problem  of  hazard,  primarily  to  the  environment,  is  nonexistent  a  few  feet  from 
a  gamma  emitter:  However,  the  problem  should  still  be  considered  when  neutron  sources  are 
used  to  induce  activity  in  seawater.  This  problem  can  be  discussed  in  terms  of  the  AEC  require¬ 
ment  for  disposal  of  radioisotopes  into  sewers.  For  -^Na  and  J^C1.  the  maximum  permissible 
concentration  in  water  for  continuous  exposure  to  the  population  at  large  is  3  x  1 0~  ’  Ci  1 
(ref.  4).  This  is  50  times  the  concentration  calculated  for  the  initial  activation  of  seawater 
flowing  at  1  knot  past  a  10'  '-n/sec  source.  If  the  flow  is  so  slight  that  the  activity  is  essen¬ 
tially  distributed  by  diffusion,  the  concentration  of  -^Na  becomes  the  principal  hazard.  At 
equilibrium,  where  the  decay  rate  equals  the  rate  of  formation,  the  total  -^Na  activity  in  sea¬ 
water  from  a  1  O'  '-n/sec  source  is  0.03  Ci.  If  this  is  distributed  or  averaeed  over  3  volume  cf 
10-m-"*  or  more,  this  concentration  will  be  equal  to  or  below  the  AEC  maximum  permissible 
concentration.  The  equilibrium  period  for  -^Na  is  2 1 .6  hr:  thus,  if  the  flow  is  such  that  the 
water  is  exchanged  within  14  ft  of  the  source  at  least  once  in  21.6  hr.  the  volume  of  water  with 
a  concentration  above  the  AEC  limit  will  be  even  smaller.  Diffusion  rates  near  the  bottom 
generally  exceed  this  flow  rate  (ref.  2).  and  it  is  unlikely  that  a  neutron  source  implanted  in 
the  ocean  will  ever  constitute  a  hazard  to  the  environment  beyond  a  few  feel  of  the  source. 

BEACON  SYSTEM 

A  radioactive  beacon  system  consists  of  a  source  implanted  above  the  bottom  and  a 
detector  aboard  the  search  vessel.  The  source,  a  radioisotope  in  a  corrosion-resistant  housing, 
will  not  require  auxiliary1  power,  and  it  can  be  selected  to  last  “indefinitely”  compared  with 
probable  periods  of  use.  When  a  bare  source  is  handled  in  the  water,  immediate  contact  must 
be  avoided,  and  it  is  probable  that  a  small  shield  will  be  necessary  to  facilitate  implantment 
and  removal.  Some  form  of  shield  must  be  provided  on  the  surface  to  protect  personnel  dur¬ 
ing  transport  of  relatively  large  sources.  While  the  implantment  of  sources  in  the  bottom  is 
possible,  the  background  radioactivity  of  the  bottom  is  generally  10  to  1000  times  that  of 
open  seawater,  and  'Ire  detection  range  (above  background)  will  be  shortened  several  feet. 
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In  most  instances,  it  is  desirable  to  locate  the  source  6  ft  or  more  above  the  bottom.  This  is 
particularly  true  when  the  flow  of  seawater  containing  activation  products  is  used  to  extend 
the  detection  range. 

The  receivers  for  a  beacon  system  will  be  scintillation  detectors.  Underwater  sys¬ 
tems  using  Nal  detectors,  which  measure  either  4  in.  in  diameter  x  4-in.  thick  or  2  in.  in 
diameter  x  8-in.  long,  have  been  successfully  used  to  measure  the  natural  activity  levels  in 
seawater.  The  small  diameter  systems  arc  extremely  rugged  and  have  been  assembled  in  arrays 
capable  of  detecting  activity  changes  of  2  to  5  percent  of  the  background.  A  levei  of 
10  percent  of  the  activity  was  used  in  the  previous  range  calculations  (this  can  be 
obtained  with  a  single  detector). 

The  cost  of  a  beacon  system  using  gamma-ray  emitters  »vould  be  nominal,  less  than 
S5K  for  the  source  and  S20K  for  a  detector-ratcmeter  system  (assuming  the  R&D  costs  can  be 
spread  across  several  sy  stems).  Pricing  of  high-intensity  neutron  systems  is  uncertain.  A  pro¬ 
jected  price  for  -^^Cf  in  1980  is  S25K  for  10*  ^  n/sec  (AEC)  or  approximately  20  mg.  This 
amounted  to  the  world  supply  a  few  years  ago.  Availability  is  dependent  on  a  pending  AEC 
decision  to  produce  -^-Cf  in  relatively  iaige  quantities.  No  other  neutron  sources  appear 
satisfactory  for  this  application. 
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CHEMICAL  SYSTEMS  AS  POTENTIAL  UNDERSEA  NAVIGATIONAL  AIDS 


INTRODUCTION 

Tins  article  summarizes  a  variety  of  chemical  methods  which  might  serve  as  potential 
undersea  navigation  aids.  In  this  .yrvey.  two  general  categories  will  be  considered:  mapping 
of  the  natural  chemical  environment  and  introduction  of  chemical  markers  into  the  ocean- 
bottom  environment.  This  discussion  will  be  limited  to  those  chemical  parameters  which  can 
be  measured  in  siitt  by  remote  probes  or  samplers.  These  parameters  and  the  background  levels 
for  average  ocean  water  are  in  table  I.  It  should  be  noted  that  this  table  lists  only  those  parame¬ 
ters  for  wiiich  direct  in  siitt  measurement  systems  are  currently  available,  development  of  new 
in  <>in:  measurement  techniques  would  expand  this  list.  In  addition,  the  number  of  potential 
methods  would  he  increased  if  laboratory  -type  measurements  were  included  because  such 
methods  are  more  plentiful  and  arc  generally  more  sensitive  and  accurate. 

MAPPING  OF  THE  NATURAL  CHEMICAL  ENVIRONMENT 

We  define  chemical  mapping  as  the  mvc  of  natural  variations  in  a  given  chemical 
parameter  or  variations  in  the  distribution  of  several  chemical  parameters  as  navigational  aids. 
Either  the  seawater  or  the  bottom  sediment  can  be  mapped;  the  former,  however,  is  probably 
less  reliable  due  to  the  dynamic  nature  of  the  water  mass. 

Bottom  Seawater 

A  search  of  the  literature  revealed  essentially  no  information  regarding  horizontal 
variations  in  the  chemical  constituency  of  deep  ocean  bottom  water  over  a  relatively  small 
area.  i.e..  several  square  miles.  In  fact,  there  is  little  information  on  the  vertical  distribution 
of  chemicals  near  the  bottom.  It  is  probably  reasonable  to  assume  that  there  will  be  little  or 
no  horizontal  variation  in  either  the  major  chemical  constituents  or  the  composite  parameters 
(such  as  salinity  or  density):  however,  there  might  be  marked  variations  in  the  concentrations 
of  trace  constituents.  Unfortunately,  this  is  one  region  of  the  ocean  for  which  little  is  known 
concerning  the  chemistry. 

Bottom  Sediments 

There  is  considerable  information  in  the  literature  regarding  the  chemical  composition 
of  deep  sea  sediments.  Composition  varies  with  the  nature  of  the  sediment,  its  location,  par¬ 
ticle  size,  absorptive  properties,  etc.  Again,  however,  litt’e  Is  known  about  variability  over  a 
relatively  small  area  of  tire  bottom.  Therefore,  no  statement  can  presently  Ire  made  concerning 
the  usefulness  of  sediment  mapping  for  navigational  purposes. 
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UV  SCANS.  Because  a  number  of  compounds  that  fluoresce  under  ultraviolet  radiation 
exist  in  the  sediment,  an  ultraviolet  scan  of  the  bottom  might  be  a  useful  adjunct  to  visual 
observations. 

Summary 

Mapping  of  the  chemical  environment,  either  of  the  seawater  or  of  the  sediment,  does 
not  appear  promising.  Too  little  is  known  about  variations  in  the  chemical  parameters  over  a 
relatively  small  area  of  the  ocean  bottom. 

CHEMICAL  MARKERS 

Any  chemical  agent  that  can  be  introduced  into  the  sea  and  then  be  detected  can  poten¬ 
tially  serve  as  a  marker.  For  the  purposes  of  this  discussion,  the  markers  are  divided  into  two 
categories:  soluble  and  insoluble  chemical  systems. 

Soluble  Markers 

Soluble  markers  have  the  advantage  of  being  more  casilv  detected  than  insoluble 
markers,  but  they  are  also  subject  to  the  vagaries  of  the  transport  processes  that  oceur  in  the 
ocean  environment. 

There  are  three  possible  ways  by  which  these  markers  can  be  introduced:  instantaneous 
release:  pulsed  release  at  either  fixed  time  intervals  or  upon  interrogation,  and  continuous 
release.  Release  of  soluble  markers  will  generate  concentration  fields  whose  features  will 
depend  upon  parameters  which  describe  the  transport  processes  in  the  ocean.  (These  parame¬ 
ters  are  factors  such  as  currents,  eddy  diffusion,  and  shear  forces.)  Detailed  calculations  based 
c  n  an  ocean  diffusion  model  of  these  concentration  fields  have  been  made  bv  Mikhail  and  his 
associates  (refs.  10  and  1 1 ).  The  discussion  which  follows  is  based  entirelv  on  their  results. 

It  should  be  noted  that  the  model  used  bv  Mikhail  assumed  an  infinite  ocean  environment, 
thus,  it  is  not  strictly  applicable  to  a  situation  in  which  material  is  released  at  or  near  the  bot¬ 
tom.  Mikhail's  results,  however,  are  useful  for  this  article  because  the  vertical  dimensions  of 
the  marker  patches  are  relativelv  small  when  compared  with  the  horizontal  dimensions.  A 
second  point  to  note  is  that  the  discussion  is  based  on  calculations  made  for  ocean  conditions 
off  Cape  Kennedy  in  August  (table  2).  The  actual  conditions  in  am  given  zone  of  the  ocean 
at  any  given  time  can.  of  course,  be  quite  different:  ho wever.  the  results  of  Mikhail's  calcula¬ 
tions  are  representative  of  at  least  one  set  of  actual  conditions  encountered  in  the  ocean. 


Table  2.  Standard  Ocean  Conditions.  Based  on  Conditions  Existing 
Off  Cape  Kennedy  in  August  ( Ref.  1 0). 


Condition 

Measurement 

Eddy  diffusivitv  in  X  direction 

1400  m- 'hr 

Eddy  diffusivitv  in  Y  direction 

1400  m-/hr 

Eddy  diffusivitv  in  Z  direction 

0.47  m-/hr 

Horizontal  shear 

0.0  hr*1 

Vertical  shear 

23.S  hr*5 

Mean  current  velocity 

360  m/hr  (10  cm/sec) 

INSTANTANEOUS  RELEASE.  When  a  point  source  containing  a  quantity.  Mq.  of  a 
soluble  chemical  is  instantaneously  reieased  into  the  sea.  an  ellipsoidal  patch  is  formed  that 
expands  with  time  and  drifts  downstream  with  the  current.  The  volume  of  the  patch  within 
which  the  concentration  exceeds  the  minimum  detectable  amount.  Dq.  will  increase  to  a  maxi¬ 
mum.  Vmax.  with  time  and  then  decrease  to  zero.  Figure  1  shows  changes  in  the  horizontal 
dimensions  of  the  patch  with  time;  these  changes  were  calculated  by  Mikhail  for  standard  ocean 
conditions  (table  2)  in  which  the  ratio  of  the  amount  released  to  the  minimum  detectable  amount. 
mgd0.  was  2.2  x  109  m-L  The  shapes  in  figure  1  are  essentially  representative  of  the  overall 
shapes  because  the  vertical  dimensions  are  very  small  when  compared  with  the  horizontal  dimen¬ 
sions.  Because  the  patch  drifts  along  with  the  current,  the  center  of  each  patch  in  the  figure  will 
also  move  with  the  current  and  its  position  will  depend  upon  the  magnitude  and  the  direction 
of  the  current. 

Therefore,  although  concentration  gradients  exist  within  the  patch,  the  instantaneous 
release  method  can  be  used  for  navigation  purposes  only  if  the  magnitude  and  direction  of  the 
current  are  known  and  if  they  are  uniform  over  the  area  of  interest  during  the  patch's  lifetime. 

Thc  patch  volume  at  a  given  time  after  release  is  directly  proportional  to  the  ratio. 
Mq/Dq  The  time  of  the  maximum  patch  volume  and  the  patch’s  lifetime  are  also  related  to 
M0/D0  (but  not  linearly).  Hence,  patch  size  and  lifetime  will  be  relatively  small  if  a  smail 
amount  of  material  is  releaseo  and/or  if  the  minimum  detectable  concentration  for  the  substance 
is  high.  These  factors  are  illustrated  for  Rhodamine  B  dye  in  table  3. 


Table  3.  Instantaneous  Release  of  Rhodamine  B  Dye;  Patch  Size  and  Lifetime 
as  a  Function  of  Amount  Released  Under  Ocean  Conditions  Listed 
in  Table  2.  Minimum  Detection  Limit.  Dq.  =  0.01  ppb. 


j  Case  1 

Case  2 

M0.D0 

2.2  x  109 

! 

4.4  x  !09 

Mq 

22  kg 

1 

44  kg 

Patch  lifetime 

1~5  hr 

231  hr 

Time  of  maximum  volume 

96  hr 

1 2~  hr 

Dimensions  of  maximum  volume 

X 

21.000  m 

33.000  m 

Y 

900  m 

1.030  m 

Z 

S.3  m 

9.5  m 

From  t’n is  brief  discussion,  it  is  apparent  that  an  instantaneous  release  marker  would  he 
useful  only  in  situations  where  a  relatively  short-lived  marker  is  desired  and  where  the  water 
movement  is  both  known  and  uniform.  The  latter  requirement  can  be  difficult  to  fulfill.  Bot¬ 
tom  currents  are  complex  and  vary  with  topography,  depth,  surface  and  middepth  current 
dynamics,  tidal  forces,  and  sediment  slope  stability.  They  have  been  found  to  vary  from  0  to 
25  cm.  sec  although  they  more  commonly  range  from  a  few  tenths  to  approximately  5  cm  sec.* 
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In  addition,  the  magnitude  and  direction  of  hottom  currents  at  a  given  point  vary  with  time.* 
Thus,  the  usefulness  of  instantaneous  release  markers  appears  to  be  highly  limited. 

PULSED  RELEASE.  In  this  case,  there  is  a  series  of  insiantaneous  releases.  If  the 
time  interval  between  successive  releases  is  long  relative  to  the  patch  lifetime,  there  will  be  a 
number  of  independent  patches  which  emanate  from  the  same  point  and  behave  as  those 
described  in  the  previous  section.  If  lire  time  interval  is  short  when  compared  with  the  patch 
lifetime,  the  situation  will  be  similar  to  that  of  continuous  release.  Finally,  if  the  time  interval 
between  releases  is  neither  very  long  nor  very  short,  the  situation  will  be  confusing.  Thus, 
this  mode  of  release  will  be  unsuitable  for  Navigation  purposes. 

CONTINUOUS  RELEASE.  Wien  a  soluble  substance  is  continuously  released  into 
the  sea,  a  steady  state  will  be  established  in  which  there  are  concentration  gradients.  The 
shapes  of  the  isoconcentration  surfaces  will  depend  on  the  values  of  the  parameters  which 
define  the  transport  process.  For  a  given  ocean  state,  the  volume  of  the  patch  within  which 
the  concentration  exceeds  the  minimum  detectable  amount  is  a  function  of  the  release  rate, 
Rq.  and  the  detection  limit,  Dq.  For  the  standard  ocean  conditions  adopted  by  Mikhail  and 
icr  a  current  of  360  m/hr  ( 1 0  cm/sec).  the  patch  volume,  V,  is  giver?  by 

V  =  1.18  x  10r>  (Rq/Do)2-05  m3  . 

In  die  following  discussion,  the  ratio  Rq/Dq  is  106  m3/hr.  For  this  value  of  Rq/Dq  zt 
a  current  speed  of  360  m/hr  and  under  standard  ocean  conditions  (table  2),  the  concentration 
gradient  along  the  direction  of  the  current  at  steady  state  is  shown  in  figure  2.  Under  these 
circumstances,  the  patch  extends  more  than  3  km  downstream.  The  maximum  width  of  the 
patch  is  270  m:  the  maximum  height.  5  m:  and  the  time  io  achieve  steady  state.  2.5  hr.  This 
type  of  patcri  will,  thcrclore,  have  a  potential  range  of  3  km  and  can  be  used  to  locate  a  fixed 
point  on  the  ocean  floor.  A  current  speed  of  10  cm/sec  is  somewhat  higher  than  the  speed 
generally  encountered  near  the  bottom  at  10,000  to  20,000  ft;  currents  in  these  regions  range 
!rom  5  to  0.5  cm/scc.*  Mikhail's  calculations  show  that  the  patch  volume  increases  when 
the  current  speed  is  reduced  Fateh  shapes  as  calculated  by  Mikhail  for  three  current  velocities 
arc  shown  in  figure  3.  One  can  see  that  the  patch  dimension  in  the  direction  of  the  current 
remains  essentially  the  same,  while  the  dimension  normal  to  the  current  increases  as  current 
.speed  is  reduced.  In  fact,  the  maximum  patch  width  is  roughly  proportional  to  Y~’-.  where 
V  is  the  current  velocity. 

The  effects  of  the  release  rate.  Rq.  the  minimum  detectable  concentration.  Dq.  and  the 
current  velocity  on  the  range  of  continuous  release  markers  can  be  seen  in  figures  2  and  3.  In 
figure  2.  the  range  downstream  is  plotted  as  a  function  of  concentration  (in  units  of  Dq)  for  a 
fixed  ratio  of  Rq/Dq.  If  the  figure  is  rcplottcd  in  terms  of  Rq/Dq.  the  range  is  found  to  be 
directly  proportional  to  Rq/Dq.  Figure  3  indicates  that  a  reduction  in  current  velocity  de¬ 
creases  the  downstream  range  and  increases  the  upstream  and  lateral  ranges. 

The  amount  of  material  required,  the  release  rate,  the  cost  of  a  i  000-hr  operation,  and 
the  range  were  calculated  for  several  chemical  parameters.  The  results  are  in  Chart  3.  Section  I. 
It  is  obvious  that  a  continuous  release  marker  is  feasible  only  if  the  substance  car.  be  detected 
at  concentrations  below  0.1  ppb:  Otherwise,  a  prohibitive  quantity  of  material  will  be  required. 

DISPOSITION  OF  SOURCES.  To  cover  an  area  as  wide  as  10  mi-,  more  than  one 
marker  will  be  required.  Unfortunately,  in  the  case  of  soluble  markers,  it  is  not  possible  to 
distinguish  the  patch  of  one  marker  from  that  of  another.  An  obvious  way  to  overcome  this 


*h’  F.  Potter.  SUC.  private  communication. 
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I  cruic  Effect  of  iimcn:  Mrlocxty  on  patch  shape  fot  contmuoi.'  release  under  conditions 
listed  :n  uhle  2.  it(y  D/,  -  10*  nr’  'hr.  ( Refs.  1 0  and  1 1 ). 

disadvantage  is  to  use  a  different  chemical  in  each  marker.  For  example,  each  marker  could 
contain  a  different  dye  venose  fluorescent  radiation  wavelength  would  differ  from  those  of 
the  other  markers. 

insoluble  Markers 

T:;e  basic  advantage  of  insoluble  markers  is  that  they  are  unaffected  by  w :v.'r  motion. 
Because  chemicals  are  largely  passive,  however.  the>  cannot  be  easily  detected  h\  probes.  In 
addition,  insoluble  substances  cannot  be  easilv  collected  by  samplers  for  analysis.  Thus,  these 
markers  have  to  be  stimulated  by  some  activator,  and  the  stimulated  emission  has  to  he  de¬ 
tected  remotely.  The  efficiency  of  such  a  system  will  depend  on  the  range  of  the  activator  in 
seawater:  the  acdvator-'.o-ma.ker  geometry ■;  the  activation  efficiency:  the  ra^ge  of  stimulated 
emission  in  seawater:  and  the  marker-to-detector  geometry. 

There  are  three  possible  methods  by  which  insoluble  markers  can  be  located:  neutron 
activation,  x-ray  fluorescence,  and  photoluminescenee. 

NEUTRON  ACTIVATION.  A  neutron  source  (either  a  radioisotope  or  a  generator  I 
and  a  gamma-ray  detection  system  are  required  for  ibis  method  of  location.  In  seawater,  the 
ranee  of  neutrons  from  typical  sources  requires  that  the  source  be  pbeed  within  approximately 
!0  to  20  cm  of  the  marker.  Gamma  rays  emitted  from  neutron-activated  markers  have  energies 


of  approximately  0.5  to  1  MeV.  Half-thicknesses*  for  total  attenuation  and  total  absorption 
in  water  for  gamma  rays  in  this  energy  range  are  about  45  and  15  cm.  respectively.  Obviously, 
any  neutron-activation  system  will  have  an  extremely  short  range. 

To  obtain  information  of  navigational  value  with  neutron-activation  insoluble  marker 
systems.  extended  markers,  i.e.,  strips  laid  on  the  bottom  in  a  grid  pattern,  are  required:  the 
amount  of  material  needed  for  such  a  system  is  considerable.  For  example,  if  a  grid  covering 
a  3-km  x  3-km  square  is  formed  by  laying  strips  500  m  apart,  the  total  marker  length  is  14  km. 
Assuming  100  s  of  marker  material  per  meter,  a  total  of  14,000  kg  or  15  tona  of  material  is 
placed  on  the  ocean  floor. 

X-RAY  FLUORESCENCE.  Either  a  radioisotopic  or  electron-gun  x-ray  source  is 
needed  for  this  system.  The  emitted  radiations  arc  x-rays  of  slightly  lower  energy  than  ihose 
absorbed  by  the  marker.  Unfortunately,  the  range  of  x-rays  in  water  is  so  short  that  this 
method  is  nearly  impracticable.  For  example,  attenuation  half-thickness  in  water  for  -+0-keY 
x-rays  is  2  cm.  Thus,  the  x-ray  source  and  detector  are  placed  very  close  to  the  marker. 

PHOTOLUMINESCENCE.  Untraviolet  or  visible  radiation  can  be  used  to  stimulate 
photoluminescent  markers.  Both  UV  and  visible  light  have  a  range  of  about  30  m  in  seawater. 
Thus  either  “point"  or  extended  markers  that  photoluminesce  can  be  readily  detected  on  the 
sea  floor  by  scanning  with  UV  or  a  light  source. 


SUMMARY 


Based  on  this  survey  of  potential  chemical  methods  for  navigation  systems,  the  follow¬ 
ing  conclusions  can  be  made. 

1 .  Mapping  of  the  natural  chemical  environment  does  not  appear  promising  because 
there  is  not  enough  data  on  the  sanations  of  chemical  parameters  over  a  relative!}  small  area 
of  the  bottom. 

2.  Instantaneous  release  of  soluble  markers  has  limited  potential.  The  markers  can 
onl.  be  used  if  the  current  is  known  and  uniform  and  «!  a  short-lived  marker  is  desired. 

3.  Continuous  release  of  soluble  markers  is  a  promising  method  if  the  substance  can 
be  detected  at  concentration*  •x'.ou  0.1  ppb.  There  are  a  number  of  inexpensive  dy  e>  which 
can  be  used  for  this  pur  ^o«e:amc  :  these  dyes  are  Rhodamine  B.  Pontacy  1  Brilliant  Pink  B. 
Rhodainine  WT.  and  FU 

4.  A  fairly  wide  region  of  the  bottom  can  be  covered  by  using  several  continuous- 
release  markers,  each  of  which,  contains  a  different  substance. 

5  An  array  of  insoluble  markers,  composed  of  materials  that  photoluminesce  under 
UV  or  visible  light,  can  be  used. 

6.  There  is  little  data  on  the  chemistry  of  the  coean  bottom  at  depths  below  several 
thousand  feet.  This  information  is  essential  if  man  is  to  operate  and  work  in  this  zone  of  the 
ocean. 

It  should  be  reiterated  that  only  those  chemical  parameters  for  which  in  sim  measure¬ 
ment  systems  are  currently  available  were  considered  in  this  survey .  .V.nong  the  92  natural 
elements.  45  ore  present  in  the  ocean  at  concentrations  below  0.1  pph  and  in  si  in  measurement 
methods  are  not  available  for  r.:esi  of  them.  In  addition,  there  are  great  numbers  of  vvatei 
soluble  chemical  compounds  which  arc  nonexistent  in  seavvaicr  and.  again,  for  whiJi  in  situ 
measurement  methods  ha>e  not  been  developed.  Thus,  the  number  of  potentially  useful  chem¬ 
ical  parameters  can  be  considerably  increased  by  developing  new  anah  tical  techniques  for 
in  situ  measurements  and  by  improvement  arui.’or  modification  of  existing  techniques. 

*  Tats  n  th(  » -ro  a  wiiwm  tr  t&u  a  oc,  Jtxlf «»/  iht  rsaietstta  trsnitnt  upoi,  a  urs;  cm  *£!  r/J.V'a,  r  »  irt  thst  tntJant. 
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THE  DETECTION  RANGE  OF  EXTREMELY-LOW-FREQUENCY 
ELECTROMAGNETIC  WAVES  IN  THE  SEA 


INTRODUCTION 

The  propagation  of  electromagnetic  (EM)  radiation  in  the  sea  is  considerably  different 
from  its  propagation  in  air.  Seawater  is  an  electrical  conductor  that  severely  attenuates  TV 
waves  which  deliver  their  energy  to  the  charged  ions.  The  electric  field  of  a  plane  wave  ii. 
water  has  the  form  (ref.  1 ) 

E=  E0e'(,‘l’j)x^  ,  (1) 

where  &  is  the  attenuation  length  given  by 


when.'  fo  is  the  wave  frequency.  /*  =  4-  \  10' '  H.  m  is  the  permeability  of  seawater  and  free 
space,  and  a  =  4  (ohm  mH  is  the  conductivity  of  seawater.  In  a  distance  of  1  attenuation 
length,  the  wave's  energy  is  attenuated  by  9  dB. 

The  greatest  ranges  of  detection  will  be  obtained  for  long  attenuation  lengths  which 
occur  in  the  e x t re meiv-lo w-freq uency  (ELF)  range  ( 1  <  f  <  100  Hz). 

Tnis  article  will  predict  an  upper  bound  to  the  expected  range  for  practical  sources  m 
the  ELF  range  where  the  attenuation  length  is  25  in<  5  <  250  in.  The  range  of  detection  will 
depend  upon  the  parameters  which  aa*  outlined  below  and  examined  in  greater  detail  in  the 
following  sections. 

POWER.  Range  depends  on  the  power  radiated  b>  the  source,  which  is  assumed  to  he 
10W.  i.e..  P=  10W. 

ANTENNA  CONFIGURATION.  The  coupling  efficiency  of  the  available  power  to 
the  wave  fields  will  depend  on  the  tv  pe  and  size  of  the  antenna.  In  the  section.  Comparison 
of  Electric  Antennas,  the  common  underwater  electric  antennas  are  examined,  and  results 
show  that  the  short-circuited  coaxial  electric  dipole  antenna  couples  best  to  the  wave  fields  m 
seawater.  This  antenna  is  used  in  this  article  with  the  assumption  that  its  length  is  10  m 
Longer  antennas  have  better  coupling;  however.  10  m  is  the  longest  length  that  is  of  practical 
use.  In  the  section.  Magnetic  Dipole  Radiation  in  the  Deep  Ocean,  the  range  of  detection 
from  a  magnetic-loop  antenna  is  shown  to  be  approximately  equal  to  the  range  for  the  short- 
circuited  electric  dipole  antenna. 

FREQUENCY.  The  most  sensitive  dependence  on  frequency  is  through  the  attenuation 
length.  Noise  and  the  coupling  of  an  antenna  to  fields  in  seawater  are  also  dependent  on  fre¬ 
quency.  The  short  antennas  previously  mentioned  couple  better  to  fields  at  higher  frequencies. 


The  analysis  in  the  section,  Electric  Dipole  Radiation  in  the  Deep  Ocean ,  shows  that  the  best 
frequencies  are  in  the  ELF  range. 

NOISE.  In  the  section.  Signal  Detection  and  Noise  in  the  Sea.  the  sources  of  EM  noise 
in  the  ocean  are  examined.  It  is  concluded  that  atmospheric  noise  within  a  few  attenuation 
lengths  of  the  surface  dominates  the  results.  Below  this  depth,  detection  is  limited  by  receiver 
sensitivity.  The  ranges  obtained  are  extremely  sensitive  to  the  assumptions  concerning  the 
noise  spectrum. 

DEPTH.  At  depths  less  than  an  attenuation  length,  the  underwater  antenna  can  couple 
to  surface  waves  which  propagate  unattenuated  through  the  air.  The  main  propagation  loss  is 
attenuation  in  the  paths  between  the  source  and  the  surface  and  the  receiver  and  tre  surface. 

HEIGHT  ABOVE  BOTTOM.  Near  the  bottom,  the  submerged  antenna  couples  to  bot¬ 
tom  fields  which  propagate  with  less  attenuation  than  in  the  water.  In  the  average  bottom, 
the  attenuation  length  is  longer  by  a  factor  of  approximately  6.  The  ranges  obtained  for 
propagation  near  the  bottom  are  sensitive  to  the  bottom-conductivity  vaiue  which  can  vary 
by  an  order  of  magnitude  depending  on  tire  type  of  bottom  material. 

OBSERVATION  TIME.  The  coherent  detection  of  signals  in  noise  is  also  discussed  in 
this  article.  The  detection  of  a  signal  in  noise  is  improved  if  the  bandwidth  of  the  detector  is 
very  narrow,  which  is  equivalent  to  observing  the  signal  for  a  long  time.  It  is  assumed  that  the 
observation  time  is  a  maximum  of  100  sec.  producing  a  bandwidth  of  i0~-  Hz. 

Tlie  results  of  the  computations  are  in  tables  1  through  3.  These  ranges,  which  are 
order-of-magnitude  estimates  for  the  best  source-receiver  orientation,  represent  the  distance 
from  the  source  at  which  the  signal-to-noise  ratio  for  coherent  detection  is  unitv  for  the  obser¬ 
vation  of  1 00  sec. 

PLANE  WAVES  ABOVE  IN.  AND  BELOW  THE  SEA 

The  detection  range  of  EM  waves  from  sources  in  the  sea  is  related  to  the  characteris¬ 
tic  properties  of  plane-wave  propagation  in  the  air.  the  sea.  and  the  ocean  bottom. 

The  electric  fields  of  plane  waves  in  a  conducting  dielectric  medium  have  the  following 
dependence  on  time.  t.  and  position,  x: 

E  =  E0e^fcJt‘kx) .  (3) 

where  a  is  the  angular  frequency  and  k  is  the  complex  propagation  constant  given  K  the  dis¬ 
persion  relation  (ref.  1 ): 

k-=fcfq&;€-rjff).  (4) 

where p.  €.  and  care  the  permeability,  permittivity,  and  conductivity  of  the  medium,  respectively 

In  air.  c/we  «  1.  and  k  is  real  and  given  by 

k  =  —  =  — .  (5) 

3  Aa  c 


where  A  is  the  wavelength  in  air  and  c  is  the  speed  of  light. 
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In  and  below  the  sea  for  frequencies  below  1  MHz,<x/o;€  »  1  and  k  is  complex  and 


siven  bv 


k  =  (l  +JVSw5. 

Hie  attenuation  length.  5W  is  given  by 


Vb  y^<rw  b  *  (7) 

where  cr.%v  b  is  the  conductivity  of  the  water  or  bottom.  In  this  case,  the  electric  field  of  a  plane 
wave  has  the  form 

E  =  E0c^Kx/5  Vx^>  .  ( 8) 

At  a  distance  of  1  attenuation  length,  the  field  strength  is  down  by  4.3  dB. 

Tlie  attenuation  lengths  in  seawater  and  in  the  ocean  floor  are  plotted  as  a  function  of 
frequency  in  figure  1.  Hie  conductivity  of  seawater  is  assumed  to  be  o\v  =  4  mho /m  (ref.  1 ), 
and  the  conductivity  of  the  bottom  is  assumed  to  be  crb  =  0.1  mho/m  (ref.  1 ).  Also  plotted 
for  air  is  \I2*.  For  the  ranges,  r,  of  interest 


»  r  >  5b  >  5w  - 


Using  Maxwell's  equations  for  a  plane  wave,  it  can  be  shown  that  the  magnetic  field  is 
related  to  the  electric  field  by  the  intrinsic  impedance  of  water 


l~  _  (dfl 
t?~o- 


Hiis  relationship  is  true  for  ail  fields  which  are  more  than  a  wavelength  away  from  the  source. 
COMPARISON  OF  ELECTRIC  ANTENNAS 

Fields,  which  arc  at  a  great  distance  compared  with  the  antenna  s  size,  are  the  same  as 
fields  from  a  point  dipole  antenna  with  an  equivalent  dipole  moment  that  is  characteristic 
of  the  antenna.  The  electric  fields  are  proportional  to  this  dipole  moment.  In  tills  section,  the 
equivalent  dipole  moments  for  the  common  underwater  antennas  are  compared;  results  show 
th3t  the  short-circuited  coaxial  dipole  antenna  has  the  best  coupling  to  the  ocean  medium  3nd 
generates  the  largest  dipole  moment  and  radiation  fields. 


Biconica!  Antenna 


The  biconica!  antenna  consists  of  two  coaxial  cones  with  their  apices  together  a!  the 
origin.  Moore  (ref.  2)  has  calculated  the  equivalent  dipole  moment  and  characteristic  resistance. 


V 


For  an  input  power.  P.  the  dipole  moment  for  a  biconical  antenna  is 


2  J?  1  COS  dfi 

Pb:c  =  21  n1  cos  90  =  /  "  '  —  ==• 

J  N/s-  log  cot  0q/2 


(ID 


where  6q  is  the  cone  angle.  1  is  the  length  along  the  cone’s  surface,  and  N  -  is  the 

intrinsic  impedance  cf  seawater.  For  P  =  10  W.  1  =  3  m,  Pq  =  30  deg.  f  =  10  Hz: 


30dc" 

p. .  ~  =  5S0  A  m  . 

1  bic 


(12) 


For  0=10  deg. 


10  dca 

p  =  490  A  m  . 
bic 


(13) 


The  dependence  on  0Q  is  not  critical  for  practical  values  in  the  range  10  deg  <  0-  <  50  deg. 
Coaxial  Antennas 

The  coaxial  antenna  is  an  insulated  ware  with  seawater  as  the  outer  conductor.  1 1  has  two 
practical  configurations:  short  circuited  to  the  water  at  each  end  or  open  circuited  with  the 
insulator  covering  each  end.  Moore  analyzed  these  antennas  and  obtained  the  lollowins 
equivalent  dipole  moments  when  the  insulation  thickness  was  small  when  compared  with  the 
attenuation  length  in  water. 

For  an  input  power.  P.  the  dipole  moment  for  the  short-circuited  dipole  antenna  is 
given  by  (ref.  2) 


p  .=  a /  SPL  =  1700  Am  ',4) 

1  SL  V  UfJL 

for  P  =  1 0  W.  1  =  3  m.  f  =  1 0  Hz. 

The  dipole  moment  for  the  open-circuited  dipole  antenna  is  given  by  (ref.  2) 

p  =  \l± EL=  1480  Am  <!5) 

1  oc  y  tjfz 

for  the  above  parameters. 


Comparison  of  Electric  Antennas 

Tire  ratio  of  the  equivalent  dipole  moments  for  the  short-  and  open-circuited  coaxial 
antennas  is 


(16) 
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which  is  independent  of  all  parameters.  The  short-circuited  antenna  has  slightly  better  coupling 
with  the  radiation  fields. 

The  ratio  of  the  dipole  moments  for  the  short-circuited  coaxiai  and  biconical  antennas  is 


;  i'  '■ 


«l\»ca fia 


At  higher  frequencies,  the  biconical  antenna  will  have  a  greater  dipole  moment  and  a  larger 
radiation  field:  however,  the  losses  due  to  attenuation  at  higher  frequencies  more  than  offset 
the  increase  in  the  dipole  moment. 

These  results  show  that  the  short-circuited  coaxial  antenna  provides  the  best  coupling 
to  the  ocean  medium.  It  is  used  in  all  calculations  of  electric  dipole  fields. 

SIGNAL  DETECTION  AND  NOISE  IN  THE  SEA 

Electromagnetic  noise  in  the  sea  comes  from  three  known  sources. 

1 .  The  propagation  and  attenuation  of  natural  atmospheric  noise  at  the  depth  of 
the  receiver. 

2.  Noise  generated  by  turbulence  coupled  with  the  magnetic  field  of  the  earth. 

3.  Thermal  or  Johnson  noise  in  the  receiver  which  limits  the  receiver’s  sensitivity. 
These  sources  arc  the  limiting  factors  which  determine  the  detection  range  of  a  given 

source  configuration.  In  this  section,  the  detection  of  signals  in  a  noisy  environment  is 
examined,  and  the  sources  of  noise  and  the  depths  in  the  ocean  where  each  is  dominant  are 
discussed. 

Detection  in  a  Noisv  Environment 


Assume  that  the  signal  is  at  frequency  A  coherent  phase-locked  detection  system 
will  be  used,  and  the  signal  will  be  observed  for  a  time,  T.  After  this  time,  it  will  be  decided 
if  the  signal  has  been  detected  by  determining  if  the  detector  registers  a  signal  greater  than  the 
noise  levci,  i.e.,  if  the  signal-to-noise  level  is  greater  than  one. 

The  signal-to-noise  ratio  for  coherent  detection  is  given  by  (ref.  3) 


j2  =  |  S2M  dcu 

•L  N(Cd)  2  a 


where  N(ts 0  is  the  power  spectrum  of  the  electric  noise  in  V~/m~  Hz  and  S(cj)  is  the  field 
spectrum  of  the  signal.  Because  the  cw  signal  is  observed  for  a  time  T.  the  received  electric 
field  is  given  by 


1  -jwnt 

E(t)  =  jr:0(r.c,^)eJ  0 

(  0 


0<  t<  T 
otherwise 
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where  r  is  the  distance  from  the  source.  S(ce)  is  given  by  the  Fourier  transform  of  E(t): 


S(cj)  =  E0(r^0) 


(20) 


Tiius. 


sin~(cu-wn)T/2 


(cu-t^)- 


(21) 


Because  S~{.  t  is  peaked  at  u  =  oeg  and  is  of  width  1/T.  it  is  assumed  that  N(<u)  =  N(&>q)  is  con¬ 
stant  over  this  range.  For  T  >  1  /2=6j^,  this  assumption  is  very  good.  On  performing  the 
w  integration,  the  signal-to-noise  ratio  becomes 


d 


T 


TE0“(r,t70) 

N«V 


(22) 


Thevalueof  Eg~(r)for  marginal  detection  is  determined  by  setting  the  signal-to-noise  ratio  equal 
to  unity: 


A  similar  calculation  gives  the  magnetic  field  for  marginal  detection  in  magnetic  noise- 
H0-=^.  .24, 

wheP  M(w0)  is  the  magnetic-noise  spectrum. 

Atmsjihcric  Noise 

Lightning  is  the  main  source  of  atmospheric  FM  noise  in  tlie  1-  to  100-Hz  range.  When 
tlies:  waves  propagate  over  the  ocean,  they  arc  strongly  affected  by  the  boundary  conditions 
at  tie  surface.  Above  a  perfect  conductor,  only  the  normal  electric  field  and  the  tangential 
magittic  field  exist.  Seawater  is  not  a  perfect  conductor,  but  noise  3t  the  surface  is  consider¬ 
ing  Afferent  from  noise  over  land. 

Soderberg  and  Finklc  (ref.  4)  have  measured  ELF  atmospheric  noise  above  and  in  the 
sea.  to  their  paper,  they  discuss  the  coupling  of  no»'se  above  the  sea  to  that  below  the  surface. 
Theyhave  measured  the  noise  spectrum  in  the  5-  to  500-Hz  range  at  a  23-m  depth.  Their  data 
sill  k  extrapolated  in  this  section  to  apply  to  various  depths.  Their  data  fits  the  approximate 
turn 


10~i6  „-2h/5 


fQ<  10  Hz 


Nh(“0)  = 


(25) 


3.2  x  10'18e'2h/8 


■TT, 


10  Hz<  Tq<  100  Hz 


0 


The  maenctic-noisc  spectrum  Mj^q),  is  related  to  Njjffg)  by  the  intrinsic  impedance 
of  water  (Eq.  10) 


Mh<V  =  z§r  Nh«*o 


) 


(26) 


System  Noise 

Tne  best  receiving  system  that  can  be  constructed  will  be  limited  in  its  sensitivity  by 
thermal  noise.  This  noise  will  be  calculated  and  used  to  obtain  an  upper  limit  on  the  detection 
range.  It  should  be  emphasized,  however,  that  this  limit  in  receiver  sensitivity  is  difficult  to 
achieve  in  practice. 

The  thermal-noise  voltage  induced  in  an  antenna  is  given  by  (ref.  5) 


tj/iP 

V-  /- ,  J  W 


v  y 


.27) 


where  0  is  the  temperature  measured  in  joules.  R  is  the  input  resistance,  and  T  is  the  measuring 
time.  The  electric  field  at  an  electric  antenna  is  given  by 


1 


where  1  is  the  length  of  the  antenna. 

With  the  use  of  Eq.  23.  the  equivalent  electric-noise  spectrum  can  be  calculated: 


N  (w  )  =  te  2=iYz  =  m 

JO  ? 


1- 


for  a  short-circuited  coaxial  electric  dipole  antenna.  R  is  given  by  (ref.  2) 


R  =  iS^+R„. 
8 


where  Ru.  is  the  resistance  of  the  wire  in  the  antenna.  The  noise  spectrum  of  the  receiver  is 
smallest  when  Rw  is  negligible  and  is  given  by 
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Nj(<u)  =  Oast. 
J  21 


(28) 


At  norma!  water  temperatures  (300  °K)  and  for  1  =  10  m. 


Nj(f0)  =  1.7  x  10‘27f0  V2/m2  Hz. 


(29) 


At  10  Hz,  this  theoretical  minimum  noise  is  10^  times  lower  than  the  noise  in  the  system  used 
by  Soderberg  and  Finkle  in  their  atmospheric-noise  measurements.  The  predictions  of  range 
using  this  limit  will  be  an  upper  bound  for  the  ranges  obtained  in  actual  measurements. 

The  magnetic-noise  spectrum  for  a  magnetic  dipole  antenna  can  also  be  calculated. 

The  voltage  (Eq.  26)  induces  a  magnetic  field  in  the  loop  given  by 


'J 

V  =  Nsa-H  , 


(30) 


where  a  is  the  radius  of  the  loop  and  N  is  the  number  of  turns.  The  equivalent  magnetic-noise 
spectrum  can  be  calculated  using  Eq.  24: 


N-2TV2 


40RN*2 


Mj(<0  =  ™-=  -TT-T 4  s  -TTTJ 

j  U  til  -LL-n-a*  tti-lLS-n* 


(o-ft-Ti-vr 


For  a  magnetic  dipole  antenna,  the  input  resistance.  R.  is  given  by  (ref.  2) 
R  =  2oj2/t2ga^N2 


(31) 


(32) 


and  the  magnetic-noise  spectrum  is  given  by 


w 


_  89cr 
3«~a 


(33) 


A  comparison  with  Eq.  26  shows  that 


ii  a 


(34) 


Tlie  factor  3«~a/161  ~  1  so  tha»  the  ratio  of  magnetic  to  electric  noise  is  approximate!*  the 
same  for  both  system  and  atmospheric  noise. 
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Water  Turbulence 


| 


i 


The  noise  from  water  turbulence  is  generated  by  that  component  of  water  flow  per¬ 
pendicular  tc  the  earth’s  magnetic  field.  The  electric  field  is  given  by 


E(t)  =  vfOBjQ  , 

where  v(t)  is  the  velocity  perpendicular  to  the  earth’s  magnetic  field,  Bg. 

This  electric  field  is  polarized  in  a  direction  perpendicular  to  Bg,  and  it  has  no  magnetic 
component  because  the  field  is  in  the  same  vicinity  as  the  source.  This  noise,  which  can  be  dis¬ 
criminated  against  with  a  magnetic  dipole  antenna,  will  net  be  considered  as  an  important 
noise  source. 

Comparison 

At  shallow  depths,  atmospheric  noise  will  dominate  and  limit  the  range  of  detection. 
For  a  wave  frequency  of  10  Hz,  thermal  noise  will  equal  atmospheric  noise  at  a  720-m  depth. 
Below  this  depth,  the  range  will  be  maximum  and  limited  by  receiver  noise.  Figure  2  is  a  plot 
of  depth  as  a  function  of  frequency  for  which  the  atmospheric  and  thermal  noise  are  equal. 

ELECTRIC  DIPOLE  RADIATION  IN  THE  DEEP  OCEAN 

Radiation  from  a  shot  i-circuited  dipole  antenna  in  seawater  of  infinite  extent  will  be 
discussed  in  this  section.  Graphs  are  presented  from  which  the  deteefion  range  of  the  signal 
can  be  determined  as  a  function  of  frequency,  source  powrer.  ante  .ia  length,  noise  level,  and 
observation  time. 

Cylindrical  geometry  in  which  the  dipole  is  aligned  along  the  z  axis  at  r  =  0  will  be 
used  (fig.  3).  The  surface  and  bottom  are  assumed  tc  be  a  few  attenuation  lengths  from  the 
source.  Because  the  largest  radiated  fields  are  at  z  =  0.  the  range  for  the  receiver  will  be  calcu¬ 
lated  at  z  =  0.  The  dominant  electric  field  is  in  the  z  direction  and  is  given  by  (ref.  1 ) 


32^ 


(35) 


where 


R  =  r/5w  >  1 


(36) 


is  the  range  measured  in  water  attenuation  lengths  and  p  is  the  equivalent  dipole  moment  of 
the  antenna.  The  magnetic  field,  in  the  9  direction,  is  related  to  the  electric  field  by  the 
intrinsic  impedance 


1EJ"  eun/i 


(37) 
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I-‘b:uic  3.  Cylindrical  aoactiy. 


For  the  short-circuited  coaxial  antenna,  p  is  given  by  Eq.  14.  Hq.  35  thus  becomes 


*E. 


2  PIC^j  p"<7  c* 


!R 


4rr~ 


R2 


(3S) 


.  *> 

where  P  is  the  input  power  and  1  is  the  length  of  the  antenna.  For  marginal  detection.  Fz  ~  is 
given  by  Eq.  23: 


4*2N(wq)  _  e'-1* 
R~ 


(39) 


Tlie  right-hand  side  of  Eq.  39  is  a  function  only  of  the  normalized  range.  It  will  he  designated 
by 


,.-2R 

F:(  R)  ~  — — — -  . 

X  • 

R- 


where  the  subscript  i  designates  this  tuirctioi.  for  the  infinite  medium  case.  Tire  left-hand  side 
of  Eq.  39  is  a  function  of  frequency,  tin.*  product  TPt.  and  the  depth  through  the  dependence 
of  the  noise.  It  is  designated  by  the  function 


GjtTPl.h.  t'0;  = 


4~  Ntop,} 

TPIag-U-O 


(40) 


Using  the  atmospheric-noise  spectrum  from  Eq.  25,  the  function  G;(TPI.  h.  f0)  is 
plotted  as  a  function  of  frequency  in  figure  4  for  TP1  =  1 04  W  sec  m  and  h‘=  0, 400,  S00  m. 

Gj  is  also  plotted  for  the  deep  ocean  where  the  system  noise  of  Eq.  29  dominates.  In  addition. 
Fj(R)  is  plotted  as  a  function  of  R  in  figure  4.  To  determine  range  at  a  given  frequency,  C .  the 
R  is  found  for  which  Gj(fQ)  =  Fj(R)  which,  in  turn,  gives  the  range  in  attenuation  lengths, 
figure  1  can  be  used  to  convert  R  to  a  range  in  meters. 

The  maximum  ranges  possible,  found  by  using  the  curve  for  Gj  (system  noise),  are 
obtained  at  depths  below  those  in  figure  2.  The  curve  for  the  range  at  the  surface  (Gj(0))  is 
included  as  an  upper  bound  on  G:.  Within  a  few  attenuation  lengths  of  the  surface,  the  range 
will  be  longer  and  the  results  from  the  section.  Diftolc  Radiation  Scar  the  Surfaie .  should  be 
used.  If  the  source  and  receiver  arc  near  the  bottom,  the  ranges  will  also  be  longer  than  those 
in  this  section,  and  the  results  from  the  section.  Dipole  Radiation  Sear  the  Bottom,  should  be 
used. 

Table  1  is  a  summary'  of  ranges  obtained  in  the  ocean  far  fr«»m  the  surface  and  the 
bottom. 


Table  1.  Range  in  tiic  Infinite  Medium  Approximation. 
TP1  =  104  W  sec  m. 


Depth  (h). 

Frequency 

R(5W) 

m 

(0.  Hz 

Range,  m 

0 

10 

12.5 

1000* 

0 

’00 

16.7 

260* 

400 

4 

13.5 

1700 

400 

10 

17 

1360 

SOO 

• 

4 

16.3 

2000 

rrstem-  [ 

>1300  | 

4 

19 

2400 

noise  ! 

>  720  1 

10 

20 

1600 

limited  ( 

>  iSO  j 

100 

22 

330 

*Srt  tcrtien.  Pipnl<  Radiation  Xtar  th<  Surface 


MAGNETIC  DIPOLE  RADIATION  IN  THE  DEEP  OCEAN 

Radiation  from  a  magnetic  dipole  is  similar  to  that  from  an  electric  dipole  with  the 
electric  and  magnetic  fields  interchanged.  The  fields,  maximum  at  /.  =  0  (see  fig.  5 ).  are  gj\en 
by  (ref.  1) 


FREQUENCY  (fQ).  Hz 


Fircic  4.  FrojosiUoa  nactios  F-  sod  G-  for  initciti  tocdieo  ippiovirutjor- 


where 


R  =  r/£w  >  I 

and  m  is  the  magnetic  dipole  moment.  The  magnetic  field  is  related  to  the  electric  field  by 
the  intrinsic  impedance 


E2  am 


The  equivalent  magnetic  moment,  in.  for  a  loop  antenna  of  radium  a.  has  been  calculated  bv 
Moore  (ref.  2) 

2  _  3  c~Pa 

m  *» 

2uJ -p.-a 

Hz  then  becomes 

IHJ2  =  Ihmsz 

V  64 

Using  the  fact  that  magnetic  noise  is  given  by 

Mcw0)  =  SiL  NH)  ’  <44) 

the  following  is  obtained  for  marginal  detection: 


^2R 

R- 


(42) 


(43) 


64  N^p) 

3  PaoJ^aT 


=  Fj(R)  =-^ 

3s" 


1 

a 


Gj(TPl.  h.  fQ) . 


(45) 


For  practical  antennas,  the  factor 


and  the  ranges  obtained  are  almost  identical  to  those  obtained  for  the  electric  dipole  antenna. 
When  a  magnetic  detector  is  used  in  the  system-noise-liiniicu  region,  the  exact  equation  ;Eq.  39) 
results.  This  is  not  surprising  when  the  symmetry  between  the  electric  and  magnetic  fields  and 
the  noise  in  far-field  ranges  of  the  sources  are  considered. 
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DIPOLE  RADIATION  NEAR  THE  SURFACE 

When  the  source  and  receiver  are  at  a  depth  much  iess  than  the  range  predicted  for 
dipole  radiation  in  the  deep  ocean,  the  range  of  propagation  is  considerably  increased.  In 
subsurface-to-subsurface  propagation,  energy  follows  an  up-over-down  path,  and  the  ware 
suffers  exponential  attenuation  in  the  up-and-down  portions  of  the  path.  In  air  where  »hc 
wavelength  is  very  long,  geometrical  losses  dominate. 

The  problem  of  electromagnetic  waves  near  the  surface  has  been  investigated  by  Moore 
(ref.  2)  and  Banos  and  Wesley  (ref.  6).  Their  results  are  summarized  by  Kraichman  (ref.  1 ). 
For  propagation  in  and  above  seawater,  the  formulas  in  fable  3.7  of  Kraichman  are  applicable. 
In  this  section,  the  geometry'  of  figure  3,  with  the  air-water  interface  at  a  height,  h,  above  the 
origin,  will  be  used.  Formulas  are  given  for  the  electric  and  magnetic  dipoles  aligned  in  the 
vertical  (z)  direction  and  in  the  horizontal  (<P  =  0)  direction.  Comparison  of  the  field  formulas 
shows  the  following  characteristics  for  distances  r«  X  air. 

Vertical  Electric  Dipole  ( VED) 


For  the  VED.  the  electric  field  in  the  ^direction.  Ep  is  larger  than  the  fields  in  the 
other  directions  by  a  factor  of  r/6w  ^  1 .  Tire  only  magnetic  field  is  in  the  £  direction,  and  it 
is  related  to  Er  by 


_  < of* 

■  tt  "*  <r 

i*Vr 


This  ratio  is  the  same  as  the  ratio  of  the  noise:  thus,  the  ranges  of  both  fields  are  equal. 
Horizontal  Electric  Dipole  (HEDj 

For  the  HED.  the  horizontal  components  of  the  electric  field  dominate  the  z  component 
by  a  factor  of  approximately 

*> 

X~air 


Tire  horizontal  components  of  the  magnetic  field  dominate  Hz  by 


The  ratio  of  the  horizontal  electric  field  to  the  horizontal  magnetic  field  is  given  by 


"I 


ill,, 


Ufl 

c r  - 


Thus,  ranges  of  both  fields  an?  equal. 


VED  and  HED  Comparison 


r  • 


The  ratio  of  the  dominant  electric  fields  in  the  VED  and  HED  cases  is  siven  by 


HED  _  V 
l^r'”  VED  Xv-J 


»  1  . 


Titus,  the  HED  has  the  larger  field  at  a  distance  r.  This  orientation  provides  the  greatest  range, 
and  it  will  be  used  in  this  section  to  calculate  the  expected  range. 

Range  will  be  calculated  using  two  assumptions:  (1 )  the  source  is  a  horizontal  dipole 
aligned  along  the  direction  9  =  0.  and  (2)  the  receiver  is  at  the  same  depth  and  is  aligned  to 
receive  the  greatest  electric  field.  The  magnitude  of  this  electric  field  has  a  very  weak  depen¬ 
dence  on  c»and  is  given  by  (ref.  1 ) 


i  p  a;n'pJa  *4h/$w  frcK-o  •>  \ 

E“  =  - — P-1 —  -  (cos  9  +  sin-  o  j  . 

8a-  r6  \  4  7 


where  R  =  r/5w  and  we  assume  that  r/Xajr  «  1.  Using  the  equivalent  dipole  moment  for  a 
short-circuited  dipole  antenna  (Eq.  14)  and  choosing  <?=  jr/2.  the  following  is  obtained: 

lep  ■  .  ,47 

-2  r6 

For  marginal  detection.  |Ej-  is  given  by  Eq.  23: 


~)  4h/o 

~~N(ccq)c  " 

->  f 

TPlt^"pt~cr 


=  _J_ 


Tlie  functions  GS(TP1.  Ii.cjq)  and  FS(R).  analogous  to  Gj  and  Fj  for  the  infinite  medium  case, 
are  siven  by 

2  .  4l»/ow 

Gs(TPi.  h.  f  .)  =  NtoQ)<;  -  <49) 

TPl^'Q-p'cr 

—  - - ft 

|-S<R»  =  I;K  . 

where  s  designates  the  functions  for  propagation  near  the  surface. 

Because  these  calculations  arc  for  radiation  near  the  surface,  the  atmospheric  ne-f-c 
will  dominate  and  Ec«.  25  will  be  used  for  N(cjq). 

FS(R)  is  plotted  as  a  function  of  R  =  r/ow  >n  figure  5.  Gs(fo>  is  also  plotted  as  a  func¬ 
tion  of  frequents  for  depths  h  =  0  anti  100  in  and  TP1  =  10't  W  m  sec.  To  find  the  range  at  a 
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given  frequency,  Tq,  the  R  at  which  F$(R)  =  Gs(fo)  is  determined.  R  is  the  range  measured  in 
units  of  8W,  which  can  be  converted  into  meters  by  using  figure  1.  Tabie  2  is  a  summary'  of 
the  ranges  obtained  near  the  surface. 


Table  2.  Summary  of  Ranges  Obtained  Near  the  Surface. 


Depth  (h),  m 

Frequency  (0,  Hz 

R 

Range,  m 

0 

4 

90 

11,250 

0 

10 

170 

13,600 

0 

100 

750 

19.000 

100 

4 

63 

7900 

100 

10 

95 

7600 

100 

100 

70 

1700 

An  analysis  of  magnetic  dipoles  aligned  vertically  and  horizontally  shows  that  the  hor¬ 
izontal  orientation  has  the  larger  fields  and  that  the  ranges  obtained  are  approximately  the 
same  as  ranges  for  the  horizontal  electric  dipole.  As  previously  stated,  the  difference  in  the 
function  G  is  a  factor  of  1 61  /(3  « -a)  >  1 . 

DIPOLE  RADIATION  NEAR  THE  BOTTOM 


The  problem  of  dipole  radiation  near  the  bottom  is  similar  to  the  problem  of  radiation 
near  the  surface.  The  main  differences  are  listed  below. 

1.  The  ware  is  exponentially  attenuated  in  the  path  through  the  bottom,  while  there  is 
no  attenuation  through  the  air. 

2.  Tlie  wavelength  in  the  bottom  is  much  longer  than  the  wavelength  through  water  so 
that  the  formulas  presented  by  Kraichman  (ref.  1 .  table  3.7)  are  applicable.  The  range  near  the 
bottom  will  be  shorter  than  the  range  near  the  surface,  but  longer  than  the  range  in  the  deep 
ocean. 

The  analyses  of  the  relative  size  of  the  field  strengths  for  vertical  and  horizontal  electric 
dipoles,  similar  to  those  previously  discussed,  have  the  following  results. 

1.  The  dominant  electric  fields  from  the  VED  and  HED  are  in  the  horizontal  plane 
parallel  to  the  bottom. 

2.  The  electric  field  from  an  HED  source  is  larger  than  the  field  from  a  VED  source 
by  the  factor  r/S^  >  1.  The  best  ranges  will  be  obtained  from  an  HED  source. 

In  the  calculation  of  range  at  the  bottom,  it  will  be  assumed  that  the  source  is  a  hori¬ 
zontal  electric  dipole  aligned  along  the  direction  <?=  0.  !t  will  also  be  assumed  that  both  the 
source  and  receiver  arc  on  the  bottom  and  that  the  receiver  is  aligned  to  receive  the  largest 
electric  field.  This  field  is  given  by 


P^O^b 


32«-<r, 


w 


(50) 
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where  Rb  =  r/Sb  >  1  is  the  range  in  units  of  the  attenuation  length  in  the  bottom.  Using  the 
dipole  moment  for  the  short-circuited  coaxial  antenna  (Eq.  14)  and  using  Eq.  23  for  the  value 
of  the  electric  field  for  marginal  detection,  the  following  is  obtained: 

GfafTPl. fg^  =  Fb(Rb,<?)  ,  pi) 

where 

Gb(TPl.  f0)  = 

TPlto-)x-crbJ> 

-2Rb  /  ? 

Fb,Rb.*)  =  -^4- 

Rb'  \ 4 

Because  these  calculations  ,>re  for  radiation  near  the  bottom,  it  will  be  assumed  that 
atmospheric  noise  is  smaller  than  system  noise,  and  Eq.  29  will  be  used  for  Nfoj).  In  figure  6. 
the  function  Fb  is  plotted  as  a  function  of  Rb  for  <?  -  0.  -/ 2.  Gb  is  also  plotted  as  a  function 
of  fQ  for  TP!  =  1 04  W  sec  m.  crw  =  4  mho/m.  and  <rb  =  0.1  mho/m.  As  in  figures  4  and  5,  the 
range  at  frequency  fQ  can  be  obtained  by  finding  the  Rb  for  which  Gb(f0)  =  Fb(Rb).  The 
range  in  units  of  the  bottom  attenuation  length.  5b.  is  Rb  =  r/5b.  The  range  in  meters  can  be 
obtained  by  using  figure  1.  Table  3  is  a  summary  of  the  ranges  obtained  at  the  bottom  for 
depths  below  those  in  figure  2  where  the  system  noise  limits  the  rangc. 


.  o 
sin-s? 


R, 
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These  results  show  the  following  characteristics. 

1 .  The  longest  ranges  are  obtained  for  propagation  near  the  surface  where  the  wave 
propagates  unattenuated  through  the  air. 

2.  In  the  deep  ocean,  the  range  is  longer  near  the  bottom  where  propagation  under¬ 
goes  le«i  attenuation.  Tins  calculation  is  sensitive  to  the  local  value  of  the  bottum  conductivity. 

3.  In  the  ocean  far  from  the  bottom,  the  greatest  ranges  are  obtained  at  low  freauencies 
and  deep  depths  where  atmospheric  noise  is  greatly  attenuated. 
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UNDERWATER  PROPAGATION  OF  LIGHT 


INTRODUCTION 

Increased  absorption  (wavelength  dependent)  and  scattering  by  water  and  its  particulate 
and  dissolved  impurities  differentiate  the  visibility  problems  encountered  underwater  from 
those  encountered  in  the  atmosphere.  Much  information  on  the  effects  of  water's  absorption 
and  scattering  on  image  degradation  is  available  (refs.  I  and  2).  Within  the  present  state-of-the- 
art.  analytic  results  are  available  where  the  distance  from  the  object  to  the  observation  point  is 
small:  however,  at  distances  which  are  large  compared  with  a  scattering  length,  one  is  often 
forced  to  use  the  Monte  Carlo  ray-tracing  methods. 

The  present  article  is  concerned  only  with  determining  the  maximum  detection  range 
of  an  underwater  light  source.  This  maximum  detection  range  will  place  an  upper  bound  on 
the  distance  obtainable  in  underwater  viewing  and  the  maximum  operating  distance  of  a  light- 
based  detection,  navigation,  or  communication  system. 

In  the  following  portions  of  this  article,  we  will  consider  various  light  sources  and  show 
that  lasers  appear  to  be  most  efficient  for  producing  maximum  on-axis  illumination.  The 
detection  of  light  will  be  shown  to  be  limited  by  the  requirement  that  approximately  10  pho¬ 
tons  reach  our  detector  in  the  required  time  interval.  Calculations  will  then  be  made  assuming 
first  distilled  water  and  then  clearest  ocean  water  for  the  case  of  no  background  illumination. 
Finally,  it  will  be  shown  that,  for  depths  less  than  approximately  4000  m.  both  hioiuminescencc 
and  daylight  filtering  in  from  the  surface  further  limit  the  maximum  attainable  range. 

SOURCES 

Three  types  of  sources,  which  appear  most  applicable  to  the  underwalei  environment, 
will  be  considered.  Tne  first  two  are  lasers  which  have  their  power  coiKentiatcd  in  wavelengths 
within  the  region  of  water's  greatest  optical  transmission,  i.e..  between  450C  and  5500  A.  The 
characteristics  of  these  lasers,  the  neodymium-doped  frequency -doubled  y  ttrium-aluminum- 
gamet  (Nd.YAG)  laser  and  the  argon  laser,  are  listed  in  table  1.  The  characteristics  listed  for 
the  YAG  laser  are  better  than  the  state-of-the-art.  and  assume  the  near-term  solution  of  the 
collimation-conversion  problem.  Hie  most  significant  differences  between  the  two  lasers  arc 
listed  below-. 

_  _  © 

1.  The  availability  of  energy  near  4750  A.  which  is  often  the  wavelength  of  least  atten¬ 
uation  by  clearest  ocean  water,  lor  the  argon  laser. 

2.  The  large  peak  powers  available  in  the  YAG  laser. 

The  latter  difference  becomes  important  when  strong  background  flux  limits  the  range 
of  the  laser.  Both  lasers  have  a  similar  total  output:  5  J/sec  for  the  YAG  laser  and  5  0  J.  sec  fo« 
the  argon  laser.  These  lasers  also  have  the  advantage  of  a  high  degree  of  collimation.  but  the 
offsetting  disadvantage  of  low’  efficiency. 

Hie  third  source,  a  thallium-iodide  (HI)  doped  mercury'  arc  lamp,  offers  the  advantage 
of  a  high  efficiency  output  between  4500  and  5500  A,  but  it  has  a  low  inherent  collimation. 


Efficiencies  of  16  percent  have  been  reported  (ref.  3)  for  Til  doped  mercury  lights,  and  the 
bmps  have  been  constructed,  with  a  reduced  efficiency,  with  arcs  as  short  as  1  mm  (ref.  4). 

For  comparison  purposes,  the  availability  of  a  'ill  doped  mercury  lamp  with  an  efficiency  of 
1 6  percent  and  an  arc  length  of  1  mnt  will  be  assumed. 

A  comparison  between  the  on-axis  illumination  of  the  Til  doped  mercury'  lamp  and  the 
argon  laser  can  easily  be  made  for  their  in-air  performance.  At  a  distance  R  from  either  source, 
the  on-axis  irradiance,  H,  will  be  approximately  given  by 


H  =  PE 


i  i 

rrR-0- 


<n 


where  P  is  the  source’s  input  power,  E  is  the  efficiency,  and  6  is  the  divergence  half-angle.  By 
inserting  a  two-lens  collimator  in  front  of  the  laser’s  output,  0  becomes 


e  = 


A1 

a2 


(2> 


where  Aj  is  the  diameter  of  the  entrance  pupil.  A?  is  the  diameter  of  the  exit  pupil,  and  0-  is 
the  initial  divergence  angle.  For  the  Til  doped  mercury  lamp  placed  :n  front  of  a  parabolic 
reflector,  we  obtain  a  beam  divergence  half-angle  of 


0  = 


(3) 


where  1  is  the  length  of  the  bmp’s  arc  and  is  the  radius  of  the  reflector.  Inserting  Eqs.  2 
and  3  and  using  the  ratio  of  ihc  laser’s  illumination.  Hj_.  to  that  of  the  TH  source.  HT1I-  we 
have 


htii  Eni  0.~Al~ 


where  we  have  assumed  equal  power  inputs  and  equal  exit  pupils.  Inserting  the  values  = 
10-3.  E-JU  =  0.16. 1  =  1  mm,  0j  =  0.45  x  1 0*3,  and  A]  =  2.5  mm.  we  obtain  H|_  =  2  x 

io4htij. 

The  illumination  in  the  forward  direction  in  air  is  then  2  x  10*  times  greater  than  that 
produced  by  the  Til  lamp.  In  water,  however,  this  ratio  will  not  be  as  great  because  of  multiple 
scatt  -ring  effects.  These  effects  affect  the  laser  and  711  beam  differently  because  the  Til  beam 
will  have  a  larger  divergence  angle  and  a  greater  proportion  cf  the  scattered  flux  will  not  be 
scattered  out  of  the  bea..i  for  equally  sized  optics.  This  means  that  the  effective  coefficient. 
orAff.  will  be  closer  to  the  true  attenuation  coefficient  for  the  laser  than  for  the  Til  bulb.  This, 
in  turn,  implies  that  the  on-axis  illumination  will  be  attenuated  more  rapidly  for  the  laser. 
Because  of  the  lack  of  measurements  or  Monte  Carlo  d3ta  for  the  distances  in  which  we  are 
interested,  this  discussion  will  be  confined  to  the  lasers  listed  in  table  1.  It  will  be  assumed 
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that  these  lasers  are  terminated  in  a  l-:n.  collimator  so  that  the  beam  divergence  half-angle  is 
0.45  x  1 CT4  rad.  In  the  next  section,  an  argument  will  be  given  to  show  that  at  the  ranges  and 
observation  angles  of  interest,  the  laser  beam  is  attenuated  as  e*0^.  where  a  is  the  true  attenu¬ 
ation  coefficient. 

RANGE  IN  DISTILLED  WATER 

In  this  section,  we  will  calculate  the  ranee  of  the  lasers  through  distilled  water,  which 
must,  of  necessity,  be  an  upper  bound  to  ranges  obtainable  in  natural  water.  It  has  been  shown 
that  filtered  natural  water,  however,  behaves  like  distilled  *vater  (ref.  5).  in  this  discussion, 
the  data  of  Clarke  and  James  will  be  considered  as  representative  of  distilled  water's 
attenuation.  It  will  be  also  assumed  that  the  detector  is  capable  of  detecting  single  photons 
(this  will  be  discussed  :n  the  next  section).  For  the  present  calculations.  5  sec  will  be  the 
decision  time,  i.c.,  it  will  be  required  that  the  mean  photon  arrival  rate  be  0.2  photons/sre- 

A  beam  with  an  initial  circular  bea.nwklih  radius.  .  and  a  half-diveracnce  ansle.  Q. 
will  produce  in  air  a  spot  of  size  r  =  vR -6-  -r  rj-  at  a  distance  R  from  the  source.  If  Fq  is  ihe 
flux  emitted  by  the  light  source.  Ihe  on-axis  intensity  will  be  given  hv 


-(R-0"-r  rp 


Using  a  receiver  with  an  apc.Zure  of  radius  r->.  the  Iot3l  flux  collected  ;n  t  see  w.H  be 

F  =  F0r-~l 
■»  ->  *» 

R~ff~  +  r," 


The  abcr»e  formula  is  somewhat  erroneous  when  the  spot  size  is  smaller  than  the  receiver  arer- 
f.ac:  however,  with  the  similarly  sized  optics  used  in  this  article,  this  error  i;  negligible. 

Because  of  scattering  ai.d  absorption  losses,  the  flux  in  water  is  less  than  the  value 
given  ir.  Eq.  6.  and  because  of  multiple  scattering  effects,  it  has  yet  to  be  determined  whether 
the  beam  is  exponentially  attenuated  by  the  prcuuct  of  the  range  and  the  true  attenuation 
coefficient.  Duntlcv  (ref.  6i  hus  -lone  same  work  3n  this  problem.  His  figure  5  is  flatted  in 
different  form  as  figure  1  i.t  this  article.  The  ordinslc  is  ’.he  solid  angle  subtended  by  the 
receiver  divided  by  2rr.  and  the  ahcis*a  is  the  range  in  iMenuation  lengths  for  which  using 
the  true  attenuation  coefficient  gives  a  wrong  answer  by  a  factor  of  30.  The  narrower  the 
tldd-of-victt .  the  longer  the  range  for  which  the  use  of  the  true  attenuation  coefficient  is 
justified.  Using  a  straight-line  extrapolation  ^trough  the  last  two  points  of  figure  i  and  using 
tire  present  beam  divergence  solid  angle-  fl/2?  -  1 .05  s  ItT^.  we  find  that  the  use  of  the  true 
attenuation  coefficient  is  justified  to  a  range  of  approximately  34  attenuation  lengths.  For 
our  calculations,  we  expect  die  iruc  attenua*ion  coef*h.k.i.t  fa  be  valid  over  a  romewiiat  greater 
range  for  the  following  reasons. 

I .  We  used  a  linear  extrapolation  of  figure  1 .  whereas  the  actual  cum  exhibits  upward 
curvature. 


2.  The  waters  which  are  used  in  this  article  have  a  lower  scattering-to-absorption  ratio 
than  that  reported  by  Duntley. 

Because  of  these  factors,  the  total  flux  intercepted  in  water  will  be  given  by 

F0r?VoR 

F  =  ~~  - 5-  (7) 

(R202+rf) 

Throughout  the  rest  of  this  article,  t  =  5  sec,  6  =  4.5  x  10*5  rad,  and  r?  =  1.0  in. 

For  distilled  water,  the  data  of  Clarke  and  James  shows  that  for  the  argon  ion  laser 
o  =  0.019/m  at  4765  A  and  ct=  0.047/m  at  5300  A  for  the  frequency-doubled  YAG  laser. 

In  the  5-sec  observation  time,  the  photon  flux,  FqT,  is  1.2  x  10*9  photons  for  the  4765-A  line 
of  the  argon  ion  laser  and  6.68  x  10^  photons  for  the  frequency-doubled  YAG.  Inserting 
these  values  in  Eq.  7  and  requiring  that  F  =  1  photon  in  the  5-sec  observation  interval,  we  find 
that  the  range  of  the  argon  laser  is  2170  m  and  the  range  of  the  YAG  laser  is  947  m.  The  next 
section  will  show  that  the  ranges  actually  obtainable  in  the  ocean  environment  are  less  than 
those  calculated  for  distilled  water. 

DETECTORS 

Photomultiplier  tubes  are  among  the  most  sensitive  detectors  of  light  in  the  visible 
region.  One  disadvantage  of  these  detectors  is  the  *‘dark  noise,”  i.e..  the  output  current 
observed  at  the  photomultiplier  anode  in  the  absence  of  any  light  input.  In  high  quality  photo¬ 
multipliers.  the  dark  current  can  be  reduced  by  cooling.  For  example,  data  on  the  RCA  type 
1P21  shows  that  the  photocathode  dark-current  emission  is  approximately  one  photoelectron 
per  sec  at  -1 50  °C  (ref.  7).  Because  the  residual  dark  current  after  cooling  mainly  consists  of 
multiple  electron  pulses  while  the  signal  at  low  light  levels  is  mainly  single  electron  pulses,  the 
noise  can  be  additionally  reduced  by  passing  the  output  of  the  photomultiplier  through  a 
pulse-height  discriminator.  Due  to  the  small  value  of  the  dark  current  after  cooling  3nd  the 
possibility  of  discrimination  against  the  remsinir.g  dark  current,  wr  vviii  neglect  the  photomul¬ 
tiplier  dark  current  in  the  calculations. 

One  source  of  noise  which  must  be  considered  for  any  realistic  laser-range  calculation 
is  the  noise  associated  with  the  statistical  distribution  cf  the  photon  arrival  times.  For  a 
Poisson  distribution,  the  distribution  of  photoelcctrons  will  also  be  Poisson  (ref  S)  w  ith  a  mean 
emission  rate  equal  to  the  quantum  efficiency’  of  the  photocathode  times  the  photon  arrival 
rate.  In  the  portion  of  the  spectrum  in  which  we  are  interested,  a  typical  value  for  the  quantum 
efficiency  is  20  percent.  To  have  a  90  percent  probability  of  at  least  one  photoclcctron  being 
emitted  during  the  observation  interval,  a  mean  photoclcctron  emission  rale  of  2.3  photo¬ 
electrons  or  a  photon  arrival  rate  of  1 1.5  photons  is  required  during  the  observation  interval. 
Thus,  a  value  of  1 1.5  will  be  used  for  F  in  Eq.  7  in  the  calculations  in  the  next  section. 

RANGE  IN  CLEAREST  OCEAN  WATER 


In  addition  to  the  attenuation  produced  by  distilled  water,  occar  -a ter  contains  particu¬ 
late  matter  and  dissolved  organic  substances  which  contribute  to  the  t<  -altering  and  absorp¬ 
tion.  Burt  (ref.  9)  has  taken  extensive  data  in  the  tropical  Pacific  on  the  variation  of  beam 
attenuation  with  wavelength,  depth,  and  location.  For  the  present  calculations,  we  will  use 


his  data  for  the  clearest  water  he  observed.  This  water,  found  at  450  m,  has  a  transmittance 
relative  to  distilled  water  which  rises  from  98  percent  at  4500  A  to  99  percent  at  5500  A. 

Burt’s  clearest  water  is  thus  slightly  clearer  than  the  98  percent  transmission  for  low  ocean 
water  given  by  Sverdrup,  Johnson,  and  Fleming  (ref.  10).  Using  Burt’s  clearest  water,  we 
obtain  attenuation  coefficients  of  a  =  0.039  for  the  4765-A  line  of  the  arson  laser  and 
a  =  0.06  for  the  5300-A  output  of  the  YAG  laser.  Inserting  these  values  into  Eq.  7,  we  obtain 
a  range  of  1030  m  for  the  argon  laser  and  710  m  for  the  frequency-doubled  YAG  laser.  These 
ranges  were  calculated  on  the  basis  of  no  background  illumination.  The  effect  of  a  light  back¬ 
ground  is  discussed  in  the  next  section. 

RANGE  IN  THE  PRESENCE  OF  BACKGROUND  FLUX 

Ambient  background  flux  makes  the  decision  concerning  !:ie  presence  of  the  laser  more 
difficult,  thus  reducing  the  range  below  that  given  in  the  last  section.  The  sources  of  background 
illumination  are  blackbody  radiation  of  the  water  medium,  skylight  which  penetrates  to  the 
laser's  operating  depth,  and  bioluminescence. 

The  total  number  of  photons  per  unit  area,  per  unit  wavelength  per  second,  emitted  by 
a  blackbody  can  be  expressed  as  (ref.  1 1 ) 


NA  = 


A4(ce/kT  ,) 


(8) 


where  c  is  the  velocity  of  light  =  3  x  10*®  cm/sec,  A  =  5300  A  =  5.3  x  10°  cm,  €  is  the  energy 
of  a  5300-A  photon  =  2.34  eV,and  kT  =  1/40  eV  at  room  temperature.  Inserting  these  values 
in  Eq.  S,  wc  find  NA  =  5.34  x  10"^  photons/scc  cnP.  Using  a  detector  area  of  ir(l  in.)-  = 

20.3  cm ^  and  a  100  A  =  1 0^  cm  acceptance  band,  we  find  a  photon  flux  due  to  blackbody 
radiation  of  1 .08  x  10~^  photons  per  second.  This  flux  is  extremely  small:  thus,  the  photon 
flux  caused  by  blackbody  radiation  is  completely  negligible. 

Qarkc  and  Denton  (ref.  1 2)  give  a  schematic  diagram  for  ine  penetration  of  light  into 
the  sea.  Their  figures  arc  for  a  clear  sunny  day  and  for  the  clearest  ocean  water  (diffuse  attenu¬ 
ation  coefficient  =  0.033/m).  Their  data  will  be  used  assuming  a  constant  spectral  distribution 
of  flux  in  the  region  between  4500  and  5500  A  and  neglecting  orientation  effects.  In  using 
data  for  brightest  daylight,  we  will  be  considering  the  maximum  possible  limitation  on  range 
due  to  penetration  of  surface  light.  This  limitation  extends  to  about  650  m.  at  which  depth 
bioluminescence  becomes  the  dominant  background. 

Bioluminescence  is  the  major  source  of  background  light  at  all  depths  on  dark  nights 
and  at  depths  below  about  650  m  (the  latter  is  valid  even  on  the  brightest  days  in  the  dearest 
water).  The  extensive  data  in  reference  12  shows  that  illumination  due  to  bioluminescence  has 
average  values  of  10^  jiW/cm-  at  650  m  to  ab  *  10*9  ptW/cm-  at  3750  m.  Bioluminesccnt 
activity  thus  decreases  with  depth,  but  not  necess^  ly  monolonicallv.  The  maximum  bio¬ 
luminescence  observed  at  any  depth  greater  than  6:  m.  however,  does  not  appear  to  exceed 
10*5  u\V/cm2. 

For  the  detection  of  the  laser  beam  in  the  p  -'•ncc  of  background  illumination, 
coherent  detection  will  be  assumed.  This  type  of  detection  yields  the  greatest  signal-to-noise 
level,  thus  pladng  an  upper  bound  on  the  range  for  a  given  source  power  and  background 


illumination.  Coherent  detection  can  be  accomplished  by  feeding  the  output  of  the  photo¬ 
multiplier  into  a  noiseless  (i.e.,  cooled)  load  resistor  and  using  the  voltage  across  the  resistor 
as  the  input  to  the  matched  filter. 

Current  noise  in  the  photomultiplier  tube  causes  a  white  spectral  noise  density,  6(a)), 
equal  to 


d>(w)  =  prci  .  (9) 

where  ft  is  the  photomultiplier  gain,  c  is  the  electronic  charge,  and  i  is  the  time-averaged 
photocathode  current  which  is  equal  to 

7  =  Ts  +  i;  .  (10) 

where  is  and  in  arc  the  time-averaged  laser  and  background  photocathode  currents, 
respectively. 

The  instantaneous  signal  anode  current,  is.  is  given  by 

is  =  Futjc  .  (11) 

where  i)  is  the  quantum  efficiency’  (tj  =  0.2).  F  is  given  by  Eq.  7  except  that  Fq»  is  replaced 
by  Fq\  the  instantaneous  source  flux  (Fq’  =  1.5  x  10-5  photons/sec  for  the  YAG  laser  and 
2.41  x  1G  - -  photons/sec  for  the  argon  laser).  The  instantaneous  anode  current  due  to  back¬ 
ground  is 


in  =  HArpijfjLc  , 


where  II  is  the  background  illumination  in  W/cm-:  A  is  the  detector  area  i  A  =  20.3  cm-):  p  is 
the  number  of  photons/secAV  (2.41  x  10^  photon/sec/W  at  4765  X  and  2.67  x  10^  photon/ 
sec  AV  at  5300  X):  and  r  is  the  fraction  of  illumination  that  is  passed  through  any  narrowband 
filter  which  can  be  placed  in  front  of  the  photomultiplier.  Assuming  a  constant  background 
illumination  within  the  band  4500  to  5500  X.  an  ideal  20-X  filler  will  yield  r  =  0.2. 

Coherent  detection  in  the  p>esencc  of  white  noise  is  governed  by  the  equation 


A, 


d"  = 


where  s~  =  i<.~for  our  case.  T  is  the  observation  time,  and  d  is  the  signal-to*noise  ratio.  We 
choose  d  =  2.  which  will  yield  a  false  alarm  probability  of  0.2  at  a  detection  probability  of 
0.85.  Substituting  Eqs.  10.  1 1.  and  12  in  Eq.  9.  we  hare 

6(a))  =  .u-e~Tj(MF  +  HArp) .  (1- 


where  the  factor  M  is  1  for  the  argon  laser  and  1(H>  for  the  YAG  laser.  Substituting  Eqs.  1 1 
and  14  in  Eq.  13.  we  find 


TMF2t] 
MF  +  HArp 


(15) 


Setting  T  equal  to  5  sec  leaves  an  equation  for  F.  the  flux  entering  the  detector  in  terms  of 
the  background  illumination  H.  Using  the  daylight  penetration  data  of  Garke  and  Denton  and 
the  bioluminescence  data  of  Garke  and  Hubbard,  appropriate  values  for  the  ambient  back¬ 
ground  for  a  given  depth  can  be  found,  and  F  at  this  depth  can  be  determined.  The  range  at 
this  depth  is  determined  by  inserting  the  appropriate  value  of  F  into  Eq.  7  and  solving  for  R. 
a.  assumed  independent  of  depth,  is  0.039/m  for  the  argon  laser  and  0.06/m  for  the  YAG 
laser  as  in  the  zero  background  calculations  in  the  last  section. 

RESULTS  AND  DISCUSSION 

The  variation  of  maximum  laser  range  with  depth  for  bright  daylight  surface  illumina¬ 
tion  is  shown  in  figure  2  and  Chart  5.  Section  I.  At  the  deepest  depths,  the  bioluminescence 
and  daylight  penetration  are  too  weak  to  limit  the  laser's  range:  thus,  the  range  is  photon 
limited.  For  the  argon  laser,  the  range  is  photon  limited  in  depths  over  4000  m.  while  the  YAG 
laser  retains  the  photon-limited  range  to  2500  m.  The  YAG  laser  is  less  affected  primarily 
because  its  radiation  is  chopped.  3nd  the  coherent  detector  is  better  at  selecting  this  signal  from 
the  background  than  the  cw  output  of  the  argon  laser.  At  depths  between  650  and  4000  m  for 
the  argon  laser  and  at  depths  between  650  3nd  2500  for  the  YAG  laser,  the  range  is  limited  by 
bioluminesccnce.  The  650-m  upper  limit  for  bioluminescence  will  become  shallower  at  night 
when  the  flux  penetrating  the  surface  is  less.  Various  points,  plotted  on  figure  2  for  the  da>  light- 
limited  regime  between  650  and  75  m.  illustrate  the  difference  between  the  argon  and  the 
YAG  lasers,  in  (his  region,  the  range  of  the  YAG  laser  again  decreases  more  slowly  because  of 
coherent  detection.  Although  figure  2  shows  that  the  argon  laser  has  a  greater  range  at  a  75-m 
depth,  this  is  due  to  the  use  of  or  =  0.039/m.  The  possibility  of  finding  water  with  this  clarity 
near  the  surface  is  remote,  and  a  value  of  or  s  0.06/m  is  more  appropriate.  In  this  case,  the 
near-surface  daylight  range  of  the  YAG  laser  will  surpass  that  of  the  argon  ion  laser. 

The  calculations  in  this  article  show  that  laser  ranges  in  excess  of  1000  m  in  the  ocean 
environment  are  difficult  to  obtain.  To  obtain  detection  ranges  greater  than  those  previously 
discussed,  it  is  necessary'  !o  increase  the  power  input,  receiving  aperture  area,  or  observation 
time.  For  the  photon-limited  operation,  an  increase  in  any  one  of  these  by  a  factor  of  10 
increases  the  range  by  approximately  2  attenuation  lengths:  2  attenuation  lengths  are  approxi¬ 
mately  32  m  for  the  YAG  laser  and  50  m  for  the  argon  laser.  Thus,  it  appears  that  very  little 
range  vvii!  be  gained  by  changing  the  system  parameters  used  in  this  article. 

The  major  limitation  concerning  the  accuracy'  of  the  calculations  in  this  article  is 
the  uncertainty  concerning  the  correct  attenuation  coefficient  to  use  on-axis  at  long  dis¬ 
tances.  The  maximum  possible  range  of  the  effective  attenuation  coefficient  is  from  a.  the 
absorption  coefficient  of  distilled  water.  to  sr.  the  attenuation  coefficient.  If  a  is  the  effec¬ 
tive  attenuation  coefficient,  the  range  in  ocean  water  can  be  as  large  as  the  range  for 
distilled  water.  Huntley's  results  f  ref.  6).  hew  ever  show  that  attenuation  coefficients 


used  in  the  sections.  Range  in  Clearest  Ocean  Water  and  Range  in  the  Presence  of 
Background  Flux,  are  probably  more  accurate.  The  final  determination  as  to  the  correct 
attenuation  coefficient  must  come  from  some  theoretical  method,  such  as  Monte  Carlo  methods 
or  experimental  i  leasurements.  Monte  Carlo  ray  tracing  is  difficult  at  long  ranges  and  for  small 
receiver  apertures.  Duntley  (ref.  6),  however,  indicated  that  it  was  possible  for  him  to  extend 
his  small-angle  measurements  to  at  least  40  attenuation  lengths.  Based  on  the  results  of  this 
article,  such  experiments  are  desirable. 
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ACOUSTIC  TRANSPONDERS  FOR  NAVIGATION  REFERENCE  POINTS 


INTRODUCTION 

The  prominent  use  of  acoustics  in  the  underwater  environment  is  due  to  the  ability  of 
sound  to  be  propagated  over  long  distances  and,  at  the  same  time,  to  cam7  a  substantial  amount 
of  information.  Sound  is  used  in  underwater  communication,  control,  detection,  and  naviga¬ 
tion.  The  present  article  will  examine  the  range,  power,  and  frequency  tradeoffs  in  an  acoustic 
transponder  navigation  system  using  pulse-burst  signals. 

TRANSPONDER  CONSIDERATIONS 

An  acoustic  transponder  is  a  device  which  replies  with  an  acoustic  signal  when  it  is 
interrogated  with  an  acoustic  signal.  The  use  of  acoustic  transponders  in  undersea  navigation 
is  based  on  the  fact  that  the  speed  of  sound  in  water  is  nearly  constant  and  that  it  can  be  mea¬ 
sured  or  calculated  when  necessary'.  Thus,  the  time  needed  for  an  interrogation  signal  to  travel 
from  the  navigated  vehicle  to  the  transponder  and  for  the  reply  signal  to  return  to  the  vehicle 
is  a  direct  measure  of  the  range  from  the  vehicle  to  the  transponder.  By  a  suitable  choice  of 
signal  encoding,  the  interrogations  and  replies  of  several  different  transponders  can  be  dis¬ 
tinguished,  allowing  the  use  of  several  independent  reference  points  in  a  navigation  system. 

Because  acoustic  interrogation  follows  the  same  path  as  a  transponder  reply  and 
because  the  two  signals  normally  use  adjacent  frequency'  bands,  this  analysis  will  consider 
only  the  one-way  travel  path.  It  is  assumed  that  the  construction  of  both  the  transponder  and 
interrogator  are  similar,  e.g.,  receiver  sensitivity,  transmitter  power  output,  and  transducer  size. 

Transmitted  power  is  limited  by  the  power  density',  which  can  be  applied  to  the  trans¬ 
mitting  transducer,  to  no  more  than  1  W'fcm-  near  the  surface.  For  this  article,  it  has  been 
assumed  that  the  transmitting  and  receiving  transducers  are  vertical  cylinders  that  are  1  wave¬ 
length  in  diameter  and  length.  Such  a  transducer  has  a  3-dB  directional  gain  in  the  horizontal 
plane  and  a  vertical  bcamvvidth  of  1 20  deg.  The  peak  transmitted  power  is  severely  limited  at 
high  frequencies  because  the  available  area  decreases  as  the  square  of  the  frequency. 

The  power  consumption  of  a  transponder  is  in  two  modes:  listening  power,  perhaps 
30  m»V.  and  transmitting  power,  which  will  typically  be  approximately  0.2  percent  of  the  peak 
transmitted  power  because  of  the  limited  duty  cycle.  Tints,  the  total  power  consumption  of  a 
transponder  while  it  is  in  action  will  be  30  inW  to  2  W. 

The  size  of  the  transponder  is  primarily  determined  by  the  size  of  its  battery  pack 
because  the  size  of  the  electronic  components  is  much  smaller  than  the  batteries.  Over  the 
frequency'  range  considered,  the  size  of  the  transponder  might  be  1 5  cm  x  1 5  cm  x  1  m.  The 
weight  will  be  10  to  50  kg,  and  the  cost  will  be  S2000  to  S5000. 


NOISE 


Hie  fundamental  question  in  this  section  concerns  the  range  at  which  a  signal  can  be 
transmitted  so  that  it  can  be  recognized  above  the  background  noise  at  the  other  end  of  the 
path.  If  this  is  the  interrogation  signal,  the  transponder  must  recognize  its  interrogation  to 
generate  a  reply.  Similarly,  the  receiving  apparatus  at  the  vehicle  must  recognize  the  reply  from 
the  transponder  to  measure  the  acoustic  travel  time  and  thus  to  measure  the  range  to  the  trans¬ 
ponder.  The  assumption  will  be  made  that  the  signal  energy  in  the  received  band  must  be  at 
least  1 0  dB  greater  than  the  background  noise  in  the  same  band. 

In  the  general  frequency  range  of  interest  (10  kHz  to  1  MHz),  there  are  two  noise 
mechanisms.  Up  to  approximately  30  kHz,  the  dominant  noise,  decreasing  with  increasing 
frequency,  is  generated  by  wave  motion  at  the  sea  surface:  above  30  kHz,  thermal  noise, 
increasing  with  increasing  frequency,  is  dominant.  The  noise  energy  produced  in  the  receiver 
will  be  proportional  to  the  spectrum  level  of  the  external  noise  and  to  the  receiver  bandwidth. 
The  receiver  bandwidth  should  be  matched  to  the  bandwidth  of  the  signal  and  hence  to  the 
length  of  the  signal  and  the  range  resolution  of  the  system. 

CALCULATION  OF  SIGNAL  RANGE 

The  signal  energy  available  at  the  receiver  will  depend  on  the  strength  of  the  transmitted 
signal  and  on  the  loss  which  the  signal  undergoes  during  its  travel  through  the  water.  The  deci¬ 
bel  loss  of  the  signal  as  it  passes  from  the  transmitter  to  the  receiver  is  a  function  of  only  the 
medium  and  the  range;  it  is  given  by 

transmission  loss  =  20  log  R  +  aR  ,  ( 1 ) 

where  a  is  the  (frequency  dependent)  attenuation  coefficient  and  R  is  the  range.  The  allow¬ 
able  transmission  loss  is  given  by 

allowable  transmission  loss  =  71.6  +  10  log  source  pow'er 

+  directivity'  index  (source) 

+  directivity  index  (receiver)  (2) 

-  10  log  bandwidth 

-  noise  level  —  recognition  differential , 

where  71.6  is  a  unit  conversion  factor.  The  source  powrer  is  determined  by  the  transducer's 
area,  proportional  to  the  wavelength  squared,  and  the  cavitation  limit  of  1  W/cm-.  For  the 
transducers  used  in  this  article,  the  directivity  index  is  3  dB  for  both  the  source  and  the 
receiver.  A  recognition  differential  of  10  dB  is  used,  and  the  noise  levels  at  various  frequencies 
are  given  in  Chart  6,  Section  I.  For  the  present  calculations,  the  resolutions  are  1  m  at  1 0  kHz 
and  0.1  m  at  1000  kHz.  The  resolutions  used  for  the  higher  frequencies  are  considerably  lower 
than  the  maximum  resolution;  however,  in  the  present  calculation,  resolution  at  the  high  fre¬ 
quency'  is  sacrificed  for  an  increased  pulsclength  which  lowers  the  bandwidth;  this,  in  turn, 
lowers  the  background  noise  and  increases  the  range.  By  inserting  the  appropriate  values  in 
Eqs.  1  and  2  and  equating  them,  the  maximum  range  is  obtained. 


RESULTS 


Chart  6  of  Section  I  illustrates  the  variation  of  maximum  range  with  frequency.  The 
large  decrease  in  range  with  increasing  frequency  is  due  mainly  to  the  rapid  rise  in  the  acoustic 
attenuation  coefficient.  The  favorable  range-bandwidth  characteristics  are  evident,  especially 
at  frequencies  of  100  kHz  and  below. 

Figure  1  indicates  the  possible  variation  of  maximum  range  with  allowable  transmission 
loss.  A  1 0-dB  degradation  of  allowable  loss  is  possible  because  the  ambient  noise  level  is  higher 
than  the  assumed  threshold;  however,  10-dB  improvements  are  possible  by  using  directive  trans¬ 
ducers  and  signal  processing  to  increase  the  effective  transmitted  power. 
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MECHANICAL  METHODS  FOR  DEEP  SUBMERGENCE  NAVIGATION 


INTRODUCTION 

Although  mechanical  navigation  aids,  the  first  to  be  developed,  are  no!  tlx-  best  or  most 
versatile  means  of  navigating  in  the  deep  ocean,  they  are  simple,  economical,  and  reliable,  and  if 
they  are  combined  with  dead-reckoning  methods,  a  usefui  navigation  system  can  be  produced. 

This  article  will  describe  some  of  the  mechanical  navigation  aids  that  have  been  used 
aboard  deep  submergence  vehicles  such  as  the  bathyscaph  TRIESTE  (refs.  1  through  6). 

EXAMPLES  OF  MECHANICAL  AIDS 

Plastic  Markers 

In  1963.  TRIESTE  was  used  in  the  search  for  the  nuclear  submarine  USS  THRESHER, 
which  disappeared  in  approximately  8500  ft  of  water  260  mi  off  the  New  England  Coast.  At 
this  time.  TRIESTE  was  the  only  Navy  submersible  capable  of  operating  at  such  a  depth. 

One  of  the  first  navigational  difficulties  TRIESTE  encountered  was  locating  useful 
landmarks  on  the  seafloor.  A  met!  od  that  proved  fairly  satisfactory  used  144j  plastic  reference 
markers  which  were  placed  every  58  ft  as  an  1 1-line  grid.  The  markers  were  numbered  from  ! 
at  the  south  end  to  131  at  the  north,  and  each  line  was  approximate!}  750  ft  apart  in  the  east- 
west  direction  (fig.  1 ). 

The  markers,  tied  to  1 0-lb  sash  weights,  were  dropped  from  an  auxiliary  ocean  tug 
everv  6  sec.  The  tug  used  Decca  purpie  and  green  lines  to  control  the  track.  Because  neither 
the  track  nor  the  logistics  of  spacing  could  be  held  constant  at  the  surface  and  because  each 
marker  took  approximately  1  hr  to  reach  the  bottom,  the  effects  of  variable  currents  consiuer- 
ably  distorted  the  actual  landings.  These  "fortune  cookies."*  however,  were  effective  in 
identifying  the  various  debris  of  the  THRESHER,  although  2  hr  sometimes  passed  between 
sightings  of  markers. 

On  one  dive.  3fter  initiating  her  descent  with  Decca  data  supplied  by  the  support  ship. 
TRI ESTE  descended  outside  the  grid  line  and  proceeded  on  a  fixed  course  awa}  from  the  grid 
because  no  markers  were  sighted.  Thus,  although  the  grid  was  approximate!}  Ur-mi  square, 
it  was  not  large  enough  to  allow  for  error  in  portion  3l  the  surface  or  for  lateral  movement 
during  descent. 

In  *964.  a  second  set  of  plastic  markers,  which  could  alwa}s  be  read,  was  dropped  in 
the  same  area.  Both  sets  of  markers  and  the  THRESHER  debris  appeared  in  photographs 
taken  !y  cameras  rigged  on  the  TRIESTE  or  towed  by  surface  ships. 

'Thau  rario*  «3rJ  far  runt  cooties  beesztse  lies  hzd  :o  unroB  before  the  rsnsxce  could  be  rasd.  end  ojii  frfSrd 
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Reference  Lines 


Bottom  reference  lines  that  indicate  position  every  200  ft  would  be  useful  for  sub¬ 
merged  navigation.  These  lines  would  eliminate  the  possibility  of  passing  between  unsighted 
markers.  The  length  of  these  lines  should  be  at  least  two  standard  errors-of-estimate  of  surface- 
position  accuracy',  but  not  less  than  1  mi  in  the  open  ocean. 

For  lines  laid  in  an  east-west  direction,  a  vehicle  could  dive  approximately  1  mi  north 
of  the  estimated  midpoint  and  proceed  in  a  southern  direction  until  the  reference  lines  were 
sighted.  It  would  then  advance  to  a  position  indicator  on  the  line  to  locate  a  reference  point 
to  iniiiate  a  controlled  search. 

Tne  reference  points  would  ideally  be  indicated  by  aeoustic/optical  reflectors  moored 
approximately  10  ft  above  the  bottom.  The  points  could  also  be  indicated  by  pingers.  flashing 
lights,  or  transponders  that  reply  only  when  queried  with  a  proper  code  or  frequency. 

Unless  prohibited  by  the  nature  of  the  operation,  reference  lines  could  be  laid  by  a 
surface  vessel  or  aircraft.  Small  anchors  would  be  attached  to  marker  points  to  hold  the  line 
in  place  on  the  seafloor.  The  line,  however,  should  be  weak  enough  to  break,  if  fouled,  without 
disturbing  adjacent  sections. 

If  surface  ships  or  aircraft  could  not  be  used  to  lay  the  line,  it  could  be  done  by  a  sub¬ 
mersible.  A  method  would  have  to  be  devised,  however,  for  a  submersible  to  lay  a  neutrally 
buoyant  or  slightly  negatively  buoyant  line  sporting  periodic,  anchored-moored.  reference 
points.  If  each  reference  point  had  a  net  negative  buoyancy'  of  1 5  lb.  the  submersible  would 
become  positively  buoyant  by  that  amount  for  each  point  laid.  A  5000-ft  line  with  markers 
every  200  ft  would  cause  a  buoyancy1  cf  375  lb  for  the  submersible,  a  buoyancy  requiring 
compensation. 

Reference  Grooves 

Reference  grooves  can  be  made  by  furrowing  a  line  on  the  ocean  flour  with  a  specially 
designed  tool.  Markers  can  then  be  dropped  into  the  furrow  at  regular  intervals.  These  markers 
would  probably  remain  in  place  because  there  would  be  no  attached  lines  to  drag  them  out  of 
position.  They  could  also  be  lighter  than  other  markers  because  a  negative  buoyancy  of  2  lb 
would  probably  keep  them  ir«  place. 

TRIESTE  plowed  grooves  on  the  seafloor  by  dragging  a  350-lb  lead  ball  which  was 
attached  aft  by  a  short  line  for  control.  The  bail  left  a  definite  groove  wherever  it  was  dragged 
except  on  certain  rock-strewn  surfaces.  The  only  problem  the  TRIKSTE  encountered  was 
drag  os:  slopes:  at  these  instances,  the  work  proceeded  more  slowly.  Some  trails  made  during 
the  1963  operations  were  still  detectable  in  1964.  and  more  recent  photographs  from  M1ZAR 
show  that  some  of  the  bottom  marks  have  lasted  well  beyond  the  short  time  they  were  needed 
and  useJ. 

On  some  sca-fioor  terrain,  such  as  the  rock-strewn  bottom  encountered  in  the 
THRESHER  search,  scribed  or  scored  grooves  are  largely  ineffective.  Hie  trail  is  usually 
broken  or  obscured  where  gullies  or  large  rocks  appear,  and  movement  or  deposition  of  sedi¬ 
ment  near  the  shore  will  fill  the  grooves,  limiting  their  effectiveness.  This  disadvantage  can 
tv  overcome  if  nylon  or  other  lines  are  used  as  adjuncts  to  the  mark  areas  where  plowing  cannot 
be  performed. 


Wheel  Den ces 


Devices  which  trail  a  wheel  over  the  bottom  have  been  used  to  measure  the  speed  of 
advance,  which,  when  coordinated  with  gyro  headings  in  a  small  computer,  yields  a  DR  plot 
automatically.  Accuracy,  however,  is  limited  where  the  bottom  is  gullied  or  where  currents 
affect  the  vehicle  and  enter  a  component  not  sensed  by  the  wheel. 

The  Perry  Cubmarine  tested  a  swiveled,  weighted,  bicycle  wheel  rim  and  fork  that 
were  rigidly  connected  to  a  6-ft  galvanized  pipe  (ref.  7).  The  distance  traveled  was  recorded 
by  a  flexible  cable-connected  odometer.  Tire  total  cost  of  the  system  was  S19  75.  Wheels, 
however,  are  not  feasible  aids  in  rocky  terrain,  on  coral  reefs,  precipitous  slopes,  or  sticky 
bottoms. 

CONCLUSION 

Relatively  effective  mechanical  methods  are  available  for  deep  submergence  navigation, 
especially  if  they  are  augmented  by  acoustical,  optical,  and  other  attainable  systems.  Basic 
mechanical  aids,  however,  can  be  used  independently  with  competent  results.  If  submersibles 
are  designed  with  adequate  cruising  capabilities  and  means  lor  bailast-loss  compensation,  they 
can  lay  reference  lines  and  markers  without  surface  assistance,  thus  simplifying  the  operation. 
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